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A bstr act 

Wavenumber' frequency spectral analysis of a 90-day winter (Jan. 15 - 
April 14) wind field simulated by a climate experiment of the GLAS atmospheric 
circulation awdel Is made using the space-time Fourier analysis which Is 
modified with Tukey's numerical spectral analysis. Computations are also made 
to examine how the model wave disturbances In the wavenumber- frequency domain 
are maintained by nonlinear Interactions. Results are compared with obser- 
vation, especially Kao and Lee's study. 

It is found that equatorial easterlies do not show up in this climate 
experiment at 200 mb. The zonal kinetic energy and momentum transport of 
stationary waves are too small In the model's Northern Hemisphere. The wave- 
number and frequency spectra of the model are generally In good agreement with 
observation. However, some distinct features of the model's spectra are 
revealed. The wavenumber spectra of kinetic energy show that the eastward 
moving waves of low wavenumbers have stronger zonal motion while the eastward 
moving waves of Intermediate wavenumbers have larger meridional motion compared 
with observation. Furthermore, the eastward moving waves show a band of large 
spectral value in the medium frequency regime. The frequency spectra in high 
frequency regime decrease faster than observation as frequency Increases. The 
scheme proposed by Kao and Lee for the maintenance of kinetic energy spectra 
by nonlinear Interactions in middle latitudes is not applicable over the whole 
model globe because of the disappearance of equatorial easterlies. The main- 
tenance of momentum flux spectra by nonlinear Interactions In NH middle lati- 
tudes Is similar to that of kinetic energy spectra. The primary nonlinear 
Interactions of kinetic energy and momentum flux are contributed by those 




between mean zonal flow and long and medium waves with low and medium fre 
qucncies. The stationary waves do not play a significant role in the non 
linear Interactions as found in observation. 


1. Introduction 


The development of the general circulation model (CtCM) in the past two de- 
cades has reached the step where the GCN is not only used to simulate the atmo- 
spheric circulation, but also employed to make medium range forecasts and 
long-term climate study. The simulation of the GCM is conventionally verified 
against the long-term averaged spatial distribution of observed atmospheric 
circulation. In fact, numerous observational studies demonstrate that the 
life cycle of the atmospheric disturbances vary from about a week for cyclone 
waves to several weeks for planetary waves. Besides the evaluation of the GCM 
simulation in the space domain, it seems logical to assess the GCM simulation 
in the time domain, too. 

In order to accomplish the second goal, Hayashi and Colder (1977) and 
Hayashi (1974) applied the lag-correlation method to make the space-time 
spectral analysis of the wave disturbances in midlatitudes and tropics of 
the CCM of the Geophysical Fluid Dynamics Laboratory (GFDL) . The spectra] 
analysis with a lag-correJatlon method is also employed by Pratt (1979) to 
analyze the midlatitude disturbances of the GCM of the National Center for 
Atmospheric Research (NCAR) . Tsay (1974b), then, uses the space-time Fourier 
analysis proposed by Kao (1968) to examine the tropical disturbances of the 
NCAR GCM. 

The technique of the aforementioned spectral analysis in space and time 
can be classified into two methods: space-time spectral analysis (Hayashi, 

1971) and space-time Fourier analysis (Kao, 1968). The former uses a lag- 
correlation method, while the latter uses a direct Fourier transform method. 
Tsay (1974a) and Pratt (1976) show these two methods are equivalent at 
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a discrete frequency, if a frequency smoothing is applied. Recently, Kao and 
Lee (1977) revised Kao's original scheme by Introducing Tukey's (1967) spectral 
analysis technique to make extensive spectral analysis of atmospheric 
disturbances. 

The GCM developed at the Goddard Institute for Space Studies (GISS) 

(Somroerville et al., 1974) has been diagnostically analyzed in various physical 
aspects, e.g. spectral energetics (Tenenbauro, 1976). However, no attempt has 
been made to investigate the wave disturbances of the current GCM at GLAS in 
the time doinaln. In order to promote our understanding of the wave disturbances 
in the GLAS atmospheric circulation model, it is of Interest to analyze the 
wave motions of this model using modified Kao's scheme and to verify the mode»l 
simulation against Kao and Lee's (hereafter referred to as KL) results. 

Tlie detailed structure of the GLAS atmospheric circulation model is 
described in detail by Somerville et al. (1974). It is a 9- level primitive 
equation model in o-coordlnate system. The grid system is 4" latitude x 5® 
longitude. The model contains a hydrological cycle, parameterization of 
subgrid- scale cumulus convection and orography. The radiation computation 
includes the calculated cloud and water vapor distribution. Arakawa's (1966) 
scheme is applied to treat the horizontal advections. Recently, the 4th 
('rder finite difference scheme is Introduced to the GCM at Goddard Laboratory 
for Atmospheric Sciences (Kalnay- Rivas, 1977). However, the climate run 
analyzed in the present study uses the second order model. j 

Data analyzed in this study are condensed from a 4-month climate experi- 
ment (D122) using the NMC analysis of 1 January 1975 as initial field. The | 

I 

i 

! 

I 
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history output data for overy 12 hr over a period of 90 days starting from 
OOGMT, 15 January 1975, arc used. The data of the first 15 days In this climate 
experiment are not used to avoid the transition period. 

The current study consists of two parts: space-time spectral analysis 

and the maintenance of large scale moving waves by non-linear interactions. 

The discussion of these two parts is presented respectively in Sections 2 
and 3 . 


2. S pectral Atialysis 
a . Meth od 

The analysis of wavenumber-frequency spectra in this study follows KL's 
procedure in which they apply the numerical spectral analysis of Tukey (1967) 
to Kao's original space-time Fourier analysis (Kao, 1968). The detailed 
procedure of computations can be found in KL. However, a brief summary of the 
modified Kao's scheme should sufficiently serve the purpose of illustration. 
Suppose an atmospheric variable, q(X,t), be a real single-valued function which 
is piecewise continuous in an interval (0,2ti) for both longitudinal X and time 
t. The complex Fourier coefficients of this atmospheric variable, Q(k,+n), can 
be evaluated by, 


Q(k, ;^) 




i(kX+nt) . , 
q(X,t)e - 'dsdt, 

'o 


( 1 ) 


where k and n are wavenumber and frequency respectively. represents the 
eastward moving waves, while fn is westward. 

The double Fourier transform is carried out in two steps. It is evaluated 


first in the space domain. 


The* wavununber- frequency Fourier coefficient* are calculated by perforaing tlu* 
Fourier tranaform of wavenumber Fourier coef ficienta, Q^(k,t), 

Q(k,+n) - ^ p” Q^(k,t)e^^"‘dt (3) 

Tukey'a numerical spectral analysis is applied after Q^(k,t) is obtained 
from (2). Let us express 

Q(k,;^) - Q^(k,:^) + i Qj(k,;^) , 

where and are the real and Imaginary parts of Q respectively. The 
wavenumber-frequency spectrum can be calculated by 

^ Pi(k,+n)Qj^(k,+n) 1 , (4) 

where P(k,ji») ■ Py(k,+n) + 1 P^(k,j?i) is the Fourier coefficient of a variable 

p(A,t). (k,+n) is a power spectrum If p(X,t) ■ q(\,t), while E (k,+n) is 

pq — pq 

the cospectrum of p and q if p(A,t) ^ q(X,t). 

I* • Latitudinal distributions of mean zonal vel o city, kinetic energy of 
stationary and transient waves, and meridional transport of westerly 
momentum 

The Linear theory of wave motions shows that the behavior of waves criti- 
cally depends upon the zonal mean state. In order to obtain some background 
of the zonal mean flow and wave motions of the CLAS model at 200 mb, the 

latitudinal distributions of zonal mean velocity, kinetic energy and meridional 
transport of westerly momentum of stationary (n-o) and transient waves (n|<0) 
waves at this level are shown in Fig. la. The observations extracted from 
various studies for the period from December to February are shown in Fig. lb 
to verify the simulation of the GLAS model. 


The zonal mean wind of this climate run in close to the January simulation 
of Somerville et ai. (1974) and the observattoni excei^ the equatorial easterlies 
surprisingly disappear. In factu the examination of 500 mb zonal mean wind in 
this climatological experiment does not show the disappearance of easterlies 
In the model tropics. 

The model kinetic energies of the transient waves of zonal and meridional 

velocities, and are very much larger than their counterparts of 

stationary wrves, and K . The observational K peaks at the equator and 

su sv su ' ^ 

at 35*N and iuis a larger value than the mcdel's. It is obvious that the 

stationary waves are not well simulated in the model. The comparison of the 

kinetic energies of stationary and transient waves between the GLAS model and the 

observation at the 200 '.evel shows that the smaller eddy kinetic energy of 

the model is mainly due to the zoiial motion of stationary waves. 

Both K. and have maximum values occurring at the same locations of 
tu tv 

the maximum zonal mean flow in the Southern and Northern Hemisphere (hereafter 

referred as SH and NH) . The s.ime situation occurs in the observation. The model 

maximum values of in the SH is much smaller than observational values. The 

observational kinetic energy of the transient waves is almost equally distributed 

In K and in the NH (KL, 1977). This equal partition of K. and does 
tu tv ^ ' tu tv 

show in the model to some extent. It should be also noted that the model 

tu 

and are larger than observational values in the tropics. 

The latitudinal distribution of meridional transport of westerly momentum 
displayed in Fig. 1 shows clearly that the model transport is equatorward in the 
regions north of 40 °N and south of 40*S, and is poleward in the global belt 
between 40**N and 40°S. The meridional transport of westerly momentum in the 
model is mainly contributed by the transient waves. The contribution from the 
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stiitlonary wavi's la concerted constructively with the transient waves. The 
meridional transports of westerly momentum by stationary and transient waves 
In the observational study are about the same in the KMl, while the later is 
dominant In the SH. The comparison of the westerly momentum transport between 
the C'.LAS model and the observational study shows that the transport is too 
strong by model transient waves and too weak by model stationary waves in the 
Nil. In addition, the momentum transport by the model transient waves in the SH 
Is slightly weaker than observational values. 

The January (Somerville et al., I91k) and July (Stone et al., 1977) 
simulations of the GISS model have shown that the model total eddy transport of 
momentum is too weak. Although the maximum value of meridional transport of 
westerly momentum in this climatological experiment Is comparative to KL’s 
observational value, it may be inferred that the weak eddy transport of momentum 
In the January and July simulation of the GISS model is due to the small contri- 
bution of stationary waves shown in Fig. la. 

It has long been well known that eddy kinetic energy is too small and 
eddy transport of momentum is Coo weak in the numerical general circulation 
models. Manabe et al. (1970) and Wellck et al. (1971) demonstrate that the 
Increasing of the horizontal resolution of numerical general circulation can 
enhance the magnitude of eddy energies and energy conversions. This suggests 
that the Increasing of horizontal resolution may well be one of the very likely 
ways to improve the simulation of stationary waves in the GLAS atmospheric 
circulation model. In fact, it has also long been realized that stationary 
waves are generally determined by the forcing of topography and stationary heat 
sources (e.g., Derome and Wiin-Nlelsen, 1971). Tlie sophisticated general circu- 
lation models usually have realistic orography, even chough it is smoothed. 
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Tlie IncreaHe of horixoncal resolution, of course, can .Mke the topography nore 
realistic. The stationary heat sources must depend U|.on the treatment of 
diabatlc heating. In other words, the refinement of model dlabatlc heating 
scheme could be another important factor to improve the simulation of stationary 
waves . 

c . Wavenumber and frequency spectra of zonal and sicrldional kinetic energy 

In order to shed light on the wave characteristic of the GLAS mode] In 
space, the wavenumber spectral analysis is made. The wavenumber kinetic energy 
spectra of zonal and meridional motions are displayed in Fig. 2. 

The zonal kinetic energy spectra show that the eastward moving waves have 
larger energy than the westward moving waves. However, the energies of these two 
classes of waves become: closer, especially In the high wavenumber regime, as 
the equator is approached. The contrast of zonal kinetic energies between the 
eastward and westward waves is consistent with KL's study. The close comparison 
between the GLAS model and KL's observational result shows that the eastward 
moving waves of the model's low wavenumber regime, which possesses most of the 
wave energy, have larger energy content than KL values. In contrast, the 
model's westward moving waves have less energy than Che KL values. It is 
interesting to note that the zonal kinetic energy of both eastward and westward 
waves at the equator Is larger than that of KL. 

Both Che observational and model kinetic energy of zonal motion in the SH 
shown In Section 2a are mainly due to the transient waves. Desbois (1975) also 
shows that the zonal kinetic energy spectrum at 45*S decreases monotonlcally as 
wavenumber Increases. The eastward moving waves may have higher energy content. 
Since Desbois' finding does not appear In this energy spectra at 38”S shown in 
Fig. 2a, we could infer Chat the zonal kinetic energy of model eastward moving 
wave of low wavenumber regime may be smaller in the SH middle latitudes. 






The zonal kinetic energy of the model stationary wave one becomes dominant at 
IS^N and 2"N. Krlshnamurti (1971) and Krlshnamurtl et al. (1973) show that 
wave one, associated with the east-west circulation at 200 mb over the tropics, 
possesses the most significant part of wave energy and is quasl-statlonary . 

This particular feature of wave one in the tropics is also confirmed recently 
by Chen (1980) and shown very clearly In KL's result. Compared with KL's 
result, the zonal kinetic energy of this wave in the equatorial area of the GLAS 
model Is still too weak. 

In the high wavenumber regime, Tenenbaijm (1976) has shown that the overall 
NH kinetic energy spectrum of the GISS model has a slope of -2.6 between waves 
8 and 15. Fig. 2a shows that the -3 power law is applicable in high and middle 
latitudes. However, the slope of the zonal kinetic energy spectrum in the high 
wavenumber regime becomes -2 as it approaches the equator. 

Model meridional kinetic energy spectra at various latitudes are shown 
In Fig. 2b. The most striking difference of the spectral distribution of 
meridional kinetic energy from that of zonal kinetic energy 's the former 
always shows a band of maximum values over the scales In which exist the 
maximum ban cl Inlc instability (KL, 1977; Pratt, 1977, Desbols, 1975). The 
GIJV.S nwdel, in general, reproduces fairly well the meridional kinetic energy 
spectra. Hovever, careful Inspection shows that the meridional kinetic energy 
in the low wavenumber regime Is too small and too large over the band of maximum 
values. In other words, the long waves have less north-south motion, while the 
most baroclinically unstable waves have too strong north-south motion. Tlie model 
westward moving waves have more meridional kinetic energy than the model east- 
ward moving waves at the equator. The slope of meridional energy spectra in 
the high wavenumber regime is close to -3, out the slope becomes -2 near the 
equator. Furthermore, the comparisen between the zonal and meridional kinetic 
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cMiergy spectra shows that the energy content ot eastward moving waves In tlie 
high wavenumber regime of meridional kinetic energy spectra is larger than 
those of zonal kinetic energy spectra, except at 2"N and 38°S. 

In order to compare the time characteristics of model disturbances with 
the observational study of KL, the frequency spectra of zonal and meridional 
motion of the eastward and westward motion at 200 mb of the Gl^S model arc 
shown in Fig. 3. There is no pronounced difference between the frequency 
spectra of zonal and meridional motion as the wavenumber spectra of these two 
motions. The spectra of the eastward moving waves have larger energy content 
than the westward moving waves, except some low frequency components of zonal 
motion at 58"N and intermediate frequency disturbances of meridional motion at 
2®N. The difference of frequency spectra between eastward and westward motions 
decreases as it approaches the equator. Furthermore, the spectra decrease as 
the frequency Increases. These general features of frequency spectra ar. 
similar to KL's analysis. Notice that some differences exist between KL's and 
the model results. The model fre*,,uency spectra of eastward moving waves have 
ii band of large values between 5-20 cycles (90 days) ^ in the middle and high 
latitudes of the NH and 38“S. The most distinct feature is that the eastward 
moving waves peak at 7 cycles (90 days) ^ which has a period of 12.9 days in the 
frequency spectra of zonal motion at 58“N, 38®N, 18®N, and meridional motion 
at 58*N, 38*N and 38®S. In fact, Hayashl and Colder (1977) also find a' 12-day 
wave in the zonal and meridional motion at 19(1 mb of the GFDL model. They 
Interpretate the 12-day disturbance as strong cyclones. Miller (1974) also 
shows a 14-15 day vacillation of atmospheric energy. In any event, this type 
of wave disturbance does not show in the KL observational analysis. Through 
careful comparison between the frequency spectra of KL and GLAK model, one finds 
that the low frequency (< 5 cycles/ 90 days) regime of the GLAS model has less 
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iiiergy than KL’s observational spectra. Pratt (1979) also finds that both 
the UFDL and NCAK models In middle and high latitudes have less wave activity 
than the atmosphere In the planetary-scale wave domain of 1-3 week period. 

The frequency spectra of the CIJVS modeJI In the high frequency regime obviously 
follow a power law as Kl/s analysis. The slope of the model freqt mey spectra 
varies from -2 to -3, with one -4 at f»8*S for zonal motion. In other words, 
the frequency spectra of the GLAS model decrease faster with Increasing 
frequency than the observation. 

• Wav enumber and frequency spectra of momentu m tran sport 

The general feature of the momentum transport In the wavenumber domain of 
the NH in tlie study of KL, and ilayastl and Colder (1975) is: (a) the equatorial 

transport at high latitudes by stationary and transient waves ef low wavenumber 
regime; (b) the poleward transport In middle latitudes by the planetary-scale 
stationary waves and barociicically unstable transient waves and (c) the south- 
ward transport at the equator by planetary-scale stationary waves. The SH in 
Desbols' (1975) study shows a somewhat different picture: (a) the poleward 

transport In the middle latitudes is attributed to the baroclinlcally unstable 
waves and (b) the equatorward transport at high latitudes is due to Icw-wavenumber 
waves . 

The wavenumber cospectra of momentum transport of the GLAS model at 200 
mb for various latitudes are shown In Fig. 4a. The momentum Is transported 
equatorward at 58®N by wavenumber 4 and 5 of the stationary and eastward moving 
waves, while wavenumber 3 the stationary wave transports momentum northward. 
The equatorward transport by wavenumber 2 of tne stationary and eastward moving 
waves found bv Kf, is not shown In the model. The poleward transport at 38“N 
is mainly accomplished by baroclinlcally unstable eastward moving waves and 
stationary wavcnui-'ber ) aad 5. Notice that the transport by barocllnlc 
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unstable eastward moving waves Is much larger than che atmosphere in KL's 
analysis. The pronounced transports of stationary wavenumbers 1 and 5 of 
the westward moving waves are missing in the model. The poleward transport 
at 18"N is weak and is attributed to stationary wave one and cyclone-scale 
eastward moving waves. The pronoun u eo southward transport at the equator 
by stationary wavenumbers 2 and 4 oi ' le KL analysis does not exist in the 
model, while the model stationary wavenumber 1 pn /ides northward transport. 
Comparing the momentum transport in the SH of the GLAS model at 200 mb 
with Desbols* result, we can see that the maximum transport is attributed 
to the eastward moving wave 7 rather than wavenumber 5 as the observational 
study shows. In addition, the equatorward transport of momentum by the 
low wavenumber regime in Desbois* analysis disappears in the model. This 
transport is mainly contributed by the cyclone-scale eastward moving waves. 

The frequency cospectra of momentum transport at high latitudes is mainly 
performed by eastward moving waves of low frequency. In middle latitudes, 

the momentum is mainly transported poleward by the eastward moving waves, 
especially those with a period of 5, 8, and 16 days, and by the westward 
moving waves of low frequency. The transport is not significant at the 
equator by the moving waves. 

Except dt the equator, the contribution to the model momentum transport 
by the westward moving waves is not significant over most frequencies. At 
38"N, the equatorward transport is mainly executed by the eastward moving waves, 
especially with frequency 7 and 13 cycles (90 days)~^, the eastward moving waves 
of low frequency transport the momentum poleward. In the middle latitudes and 
subtropics of the NH, the poleward transport of momentum is mainly carried 
out by the eastward moving waves with peak values at 7 and 11 cycles (90 days) ^ 
at 38®N and 7 cycles (90 days) ^ at 18*N. The southward transport at the equator 


iH im)i*e than obnervi’d and performed by the westward moving waves 

ul low frequency enpeclally I, 4 and 8 cycles (90 days) \ In the middle latitudes 
and subtropics of the SH, the poleward transport of momentum Is mainly pro- 
vided by the eastward moving waves with frequency less than 20 cycles (90 days) 

The momentum transport Is not significant at 50“S. 

in sunmuiry, the momentum transport by the stationary waves of low wave- 
number and the westward moving waves of low frequency are not simulated well 
by the GIJVS model. 

e . Wajy e number- fr equency spectra of zonal and meridional kinetic ene rgy, 
and mom ent um t rans por t 

It was mentioned in tl>e Introduction that two methods were employed to make 
spectral analysis of atmospheric motions: the lag-correlation method and the 

direct Fourier transform method. Hayashl and Colder (1977) and Pratt (1977) 
use the former method, while Kao and his colleagues (Kao and Wendell, 1970; 

Wendell. 1969; Kao et al., 1970; Kao and Kuezek, 1973) use the latter method. 

The general finding of their wavenumber-frequency spectra of kinetic energy and 
inonK-'ntum is the following. In the middle latitudes, the preferred spectral band Is 
oriented from low wavenumber,- low frequency of westward moving waves to high 
wavenumber and high frequency of eastward moving waves. The wave motion of the 
atmosphere in the dominant region of the spectral density is essentially of the 
Rossby type. In the tropics, the spectral band is oriented from low wavenumber 
and frequency to high wavenumber and low frequency. 

Figs. 5-7 display the frequency-wavenumber spectra of zonal motion, 
meridional motion, and momentum transport at 200 mb and various latitudes of 
the (fl.AS model. The common features of these spectra is a preferred band of 
spectra that extends from a region of low wavenumber and low frequency of the 
westward moving waves to a region of high wavenumber and high frequency of 
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tlie eastward oiovlitts waves In the nlddLe and high latitudes in both hemispheres. 
The preferred band of spectra in the tropics extends from the region of low 
wavenumber and frequency to a region of high wavenumber and low frequency. 
These common features of the GLAS model spectra are similar to Kao's. 

The contrast of zonal spectra in both hemispheres shows that the Northern 
Hemisphere has higher wave activities. This is expected because the climate 
run of present analysis is the later winter and earlier spring. The maximum 
values of spectra at 3B*N and S8*N occur at wavenumbers 3 and 4 with frequency 
of 7 cycles (90 days) Notice that wavenumbers 7 and 8 with frequency of 
15 cycles (90 days) ^ at 38“S has a significant power signal which does not 
appear in the Northern Hemisphere. 

The meridional spectra show that the maximum spectral values arc wave- 
number 6 with frequency of 10 cycles (90 days) ^ at 38®N and wavenumber 8 with 
frequency of 15 cycles (90 days) ^ at 38*S. In addition, the westward moving 
waves at the equator have higher power than the eastward moving waves. This 
is also consistent with the analysis of Kao and Kuezek (loc cit). 

The comparison between zonal and meridional spectra shows that the 
maximum spectral values of the former spectra appear in the lower wavenumber 
than that of the latter spectra. In addition the contrast between the spectra 
in both hemispheres shows that the maximum spectral values in the NH appear in 
the lower wavenumber, while they appear In the higher wavenumber in the SH. 

The same situation also applies to the spectra of momentum transport. 

The spectra of momentum transport show significant northward transport 
at 38"N and significant southward transport at 38*S. The maximum values 
of spectra exists at wavenumber 4 with a frequency of 7 cycles (90 days) ^ 
at 38*N and wavenumber 7 with a frequency of 15 cycles (90 days) The 
interhemispheric contrast of zonal and meridional spectra also applies to 
the spectra of momentum transport. 

1 I •H'ri i f'il Hlil it f *"**" ' - 
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} . Mfl 1 n t enancc of large-sc ale movinR w aves by nonlinear t nter actions 

a . Format ion ajnl_ coinput a 

The maintenance of the moving waves by the nonlinear interactions in the 
f!LAS model will be analyzed in terms of the kinetic energy equations of zonal 
and meridional motions and the momentum flux equation. The derivation of 
these equations in the. wavenumber -frequency domain can be found in Kao (1968). 
The classification of wavenumber- frequency domain in the practical computa- 
tions follow KL. Therefore, a brief description of the formulation and compu- 
tations here should be sufficient. 

The kinetic energy equations of zonal and meridional motions and the 
momentum flux equation in the physical domain can be expressed as: 


1 3u^^ 

2 “at" 


2 

u_ 9u 

a cosi)i ax 

(ul) 


uv 3u . tan^ 2 . , , v . « 

“T u v fu(v-v ) -f uF, , 

a 3((i a g 1 

(u2) (u3) lu4) (u5) 


(5) 


1 av ^ 

2 at 


2 

uv 2 

a co8(Ji ax a a4> a “ 

(vl) (v2) (v3) 


fv(u-u^) + uF^, 
(v4) (v5) 


( 6 ) 



u 3/v v3/v 


(vul) 


(vu2) 


tan^ 

a 


/ 2 2 , 

u(u -V ) 

(vu3) 


- f(u(u-Ug) - v(v-Vg)l + vF^ uF 2. (7) 

(vu4) (vu5) 

The notations used are conventional. Let us use capital letters to designate 
the Fourier coefficients of atmospheric variables, equations (5)-(7) in the 
wavenumber-frequency domain can be written in the form (Kao, 1968); 
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!u(k,n) 

(EUU) 


|v(k,n)l^ 

(EVV) 


X 

2n 


U(-k,-n) 


1 1 


lU(J,»)U(k-J.n-.) 

“ (Ul) 


- U(-k.-n) ^ I I lu (j,m)V(k-j.n-ni)l 
" J m ^ 

(U2) 


-U(-k.-n)i ^ I I tU(J,m)V(k-J,n-m)l 
J w 
(U3) 


-U(-k,-n) I f|V(k,n) - Vg(k,n)l 

(U4) 


i 

n 


U(-k,-n)Gj(k,n) ■ 
(U5) 


( 8 ) 


T_ 

2n 


V(-k.-n) i J I (JV(J.«)U(k-l,n-m)| 

J n 

(VI) 


+V(-k,-n) i n tV^(j.m)V(k-j,n-m)l 
i "‘(V2) 


+7(-k,-n) ^ I I (U(J,m)U(k-J.n-m)) 

J in 
(V3) 

l-V(-k,-n) I f(U(k,n) - Ug(k.n))- ^ V(-k,-n)G2(k,n) • , (9) 

(V4) (V5) 
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E (k.n) 
iiv 


.X. 

2Tia 


-I I 
cos^ n 


I I JU(k-J,n-m)(U(.i,m)V(-k,-n) + V(J »m)U(-k,-n) 1 

.1 «" 

(VUl) 


^ H V(k-J.n-m)lU^(j.n)V(-k.-n) 4- V^(J ,m)U(-k,-n) ] 


1 «" 


(VU2) 


“ tan^ I I U(J,m)(U(-k,-n)U(k-j,n-m) - V(-k,-n)V(k-J ,n-m) 1 
n M *4 

’ *" (VU3) 


i f T 

{U(-k.-n)tU(k.n) - Ug(k,n)J - V(-k,-nlV(k,n) - V^Ck.n))} 

(VU4) 


^ ^ U(-k,-n)fi 2 (k.n) “ V(-k,-n)G^(k.n) 

(VU5) 


( 10 ) 


The left-handed aide of (8) to (10) la respectively kinetic energy of 
zonal and meridional velocity, and momentum flux of moving waves with 
wavenumber k and frequency n. For convenience of future discussion, particular 
notations are designated to every term on the right-hand side of these equa- 
tions. ( )1 and ( )2 represent the longitudinal and latitudinal convergence 
of either energy or momentum flux respectively by nonlinear interactions. 

( )3 is sphericity effect and ( )4 4- ( )S is a cooibination of the ageostrophic 
effect, friction and computational errors. This combining effect can be 
determined by the residual method. 
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Every nonlinear Interaction tern on the right-hand aide of (8) to (10) 
haa a double euMmatlon which involvea wavenumber and frequency. Following 
Kao's analyala, the range of wavenumber and frequency In the computation of 
this study covers -20 < j «• 20 and -90 < m 90. Since a great number of pos- 
sible interactions Is Involved In these compulations, Kuo's claaulf Icatlou 
of the wavenumber-frequency apace into some particular domains la also adopted 
for the purpose of analysis. Tables 1 and 2 show the classifications of 
wavenumber and frequency, respectively. 

Notice that the classification shown in Tables 1 and 2 can provide 28 
domains and 406 possible Interactions. A combination of a wavenumber and 
frequency classification represent a domain, for example, <k , a.) designates 
the domain of long waves moving eastward with low frequency. The interac- 
tions are specified by placing the wavenumber-frequency domains side by side, 
e.g. (k , n )(k . -n ) . 

B HI m t. 

b . Maintenance of spectral energy of moving waves by nonlinear interaction s 

Wendell (1969) and sequential studies of Kao and collaborators have 
shown that the sum of the spectral energy of six wavenumber-frequency categories 
sho%m in Table 3 can explain most of the kinetic energy associated with moving 
waves. In order to examine how the spectral kinetic energy of model moving 
waves is maintained through the nonlinear interactions, only these six categories 
of moving waves will be analyzed using equations (8) and (9). The contributions 
of various interactions for various latitudes at 200 mb of the GLAS model are 
shown In Table 3. 

It is of Interest to point out that EUU and EVV of long waves 
shown in Table 3 possess more spectral energy in the Northern Hemisphere, 
while medium waves have more spectral energy in the Southern Hemisphere. The 
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lonRltudlnal convergences of kinetic energy flux, U1 and VI, over the 
model globe supply kinetic energy to eastward moving waves, but extract 
kinetic energy from the westward moving waves. The function of U1 and VI 
at the equator at 200 mb of the GLAS model is opposite to that of KL's obaerva- 
tional analysis at SOO mb where U1 and VI extract kinetic energy from eastward 
moving waves and supply to westward moving waves. The disappearance of the 
easterly flow in the model equatorial area shown in Section 2b may cause the 
I'cversiil of longitudinal convergence of kinetic energy flux at 200 mb of the 

(HAS model. It Is also of interest to point out that UI and VI of long 
waves arc very significant in the NH, but those of medium waves are generally 
significant in the SH. This feature is consistent with the latltudiudl 
distribution of EUU and EW. 

Tlie latitudinal convergence of kinetic energy flux, U2 and V2, is gener- 
ally smaller than Ul and VI, especially the contrast between VI and V2. 
Furthermore, U2 and V2 are not exactly of opposite sign to Ul and VI as KL 
found at the 300 mb analysis. In fact, the signs of U2 and V2 are randomly 
distributed to some extent. The sphericity effect, U3 and V3, is generally 
very small and is essentially zero near the equator. 

The residual term, (U4 4- U5) and (VA + V5), which is a combination of the 
ageostrophlc effect, friction, and computational error, is the same order of 
magnitude as Ul and VI. However, (U4 + U5) and (V4 + V5) have the opposite 
function to maintain the spectral energy of moving waves as compared with Ul 
and VI. In other words, (U4 + U5) and (V4 V5) supply kinetic energy 

to westward moving waves and extract kinetic energy from eastward moving waves 
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over Che whole model globe, except for the eastward moving waves with medium 
frequency at some latitudes# Notice that the function of (DA 4- U5) and 
(VA 4- V5) Is also opposite to KL's finding at the equator where KL found 
Chat (UA 4- U5) and (VA 4- V5) supply kinetic energy to both eastward and 
westward moving waves. 

Based upon their analysis for the linear and nonlinear Interactions in 
spectral kinetic energy equations, KL proposed a scheme f-ir the maintenance of 
spectral energy of moving waves In the westerly regime of the middle latitudes, 
sliown In Fig. 8, and in the easterly regime of the tropics. Although U2 and V2 
of the GLAS model at 200 mb do not show such a regularity as KL found In their 
analysis, KL's scheme for the maintenance of spectral energy of moving waves 
in the westerly regime Is In general applieaMe to the moving waves of the 
GLAS model at 200 mb. 

It is also of interest to point out that our further examination on the 
contributions of various Interactions to the spectral energy of moving waves 
reveals that Che primary nonlinear Interactions (not shown) Is due Co those 
involving the mean zonal flow with the long and medium waves of low and Inter- 
mediate frequency. Our analysis In Section 2b shows Chat Che model stationary 
waves are weak. KL found that the stationary waves play an Important role In 
the nonlinear interactions. However, the present study shows that the role 
played by the model statlonarv waves Is not vital. 

c. Malntenanc. I'f spe ct ral momentum transport flux uf moving w a ves b y 
nonlinear interactions 

Equation (9) is used Co examine how the spectral momentum transport flux 
of moving waves Is maintained through nonlinear interactions for various 
wavenumber-frequency domains and various latitudes. The resultant values of 
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various Interactions for the six wavenumber- frequency categories described in 
Section 3b are displayed in Table 8. 

The spectral flux of momentum in various wavenumber-frequency domains 
shown in Table 8 directs toward north at 38*N, 18*N and 58*S, but toward 
south at 58*N, 2*N, 18*S and 38*8. The exception occurs to the low-frequency 
westward moving waves of long wavelength at 38*N where these waves present 
rhe southward transport of momentum flux. 

The longitudinal convergence of momentum transport, VUl, is still gener- 
ally the dominant process, and its function on the momentum transport is more 
complicated than its counterpart of the zonal and meridional kinetic energy 
discussed in Section 3b. At 38*N (midlatitude) and 18*N (subtropics) VUl 
provides the northward transport of momentum to the eastward moving waves, but 
extract it from the westward moving wavev . However, the reversed process occurs 
at 38*S and 18*8, except (k„, 4n ) at 38*8 and (k^, +n ) at 18®8. At the 
equator, VUl is gc>nerallv small, but provides northward transport of momeniuin 
to both the westward and eastward moving waves. At high latitudes of the 
SH (58“S), VUl behaves in a similar way as its counterpart at the equator. 
However, VUl at 58‘*N acts randomly as far as the direction of momentum 
transport is concerned. 

The contribution of latitudinal convergence of momentum transport, VU2, 

Is generally smaller than VUl, and Is random. In some wavenumber -frequency 
categories at high latitudes (58"N) and low latitudes (2“N and IS^N), VU2 
becomes more significant than VUl. The effect of sphericity, VU3, Is gener- 
ally small and almost zero at 2'N. The residual term which involves the com- 
bination of the ageostropldc effect, frictional and computational errors, 
has the comparative magnitude with respect to VUl and has an opposite sign 
except the wave categories belong to intermediate frequency domain. 
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Following Chclr Hchome for malntnlnlng the spectral energy of the 
moving waves by nonlinear interaction!, ageostrophlc effect, and friction, KL 
propose a similar scheme to show how the spectral momentum flux of the 
moving wave Is maintained by various processes. la KL'h scheme at 40*N, 

VUl -f VU3 supply northward spectral momentum flux to the eastward moving 
waves, while VU2 + (VU4 + VU5) extract northward spectral momentum flux from 
eastward moving waves. As for the westward moving waves, the functions of 
VUl 4- VUl and VU2 + (VU4 + VU5) are reversed. At the equator, VUl 4 VU2 
ext' act northward spectral momentum flux from both the eastward and westward 
moving waves, while (VU4 4 VU5) supply northward spectral momentum transport 
to both the eastward and westward moving waves. The comparison between Table 8 
and KL’s result shows that KL’s scheme may not be applied to the GLAS model 
» oni|) I of el V, enpeelalU at the equator where t lie eister I i I'.s >|e not 'how in the 

GLAS model. At and 18*N of the GLAS model, we move VU3 to the group 

VU2 4 (VU4 4 VU5) and obtain the top diagram In Fig. 8. That Is, VUl supplies 
northward spectral momentum flux to the eastward moving waves, while VU2 4 
VU3 4 (VU4 4 VU5) extract northward spectral momentum flux from the eastward 
moving waves. The effect of VUl and VU2 4 VU3 4 (VU4 4 VU5) on the westward 

moving waves are opposite to that oi. the eastward moving waves. This scheme 

Is reversed at 38"S and 18'S. It Is because the eddy transport of westerly 
momentum Is south-poleward. At the equator of the GLAS model, VU2 generally 
Is larger than VUl and Is not so regular In sign. This does not show in KL's 
observational analysis. In order to compare with KL's study, we also modify 
KL's scheme at the equator by moving VU2 to (VU4 4 VU5). The maintenance of 
spectral momentum flux at the equator of the GLAS model Is opposite to KL’s 
scheme In sign. In other words, VUl supplies northward spectral momentum 
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transport to both eastward and westward i?oving waves, while VU2 + (VU4 4- VU5) 
extract northward spectral momentunt transport from both eastward and westward 
moving waves. 

It should be also pointed out that the detailed analysis of the contri- 
butions from various interactions to the spectral momentum transport (not 
shown) are those between the mean isonal flow and the long and medium waves with 
low and interk..c>dlate frequency. Again, the stationary waves do not play any 
significant role in nonlinear inter tetions . 

4 . SjtimmajQT 

The space-time Fourier spectral analysis and the equations of kinetic 
energy and momentum transport In the wavenumber -frequency domain proposed by 
Kao (1968) and modified by Kao and Lee (1977) using Tukey’s numerical spectral 
analysis are employed to analyze the 200-mb wind fields of a 4-month climate 
experiment of the GLAS atmospheric circulation model. This climate run covers 
the period from January 1 to April 30, 1975. This study analyzes a 90-day 
(January 15-April 14, 1975) wind field of this climate run. Comparison 
between the present study and observation, mainly Kao and Lee's study, leads 
to the following conclusions. 

1. The 200-mb latitudinal distributions of model mean zonal wind, zonal and 

meridional kinetic energy and momentum transport of eddies agree fairly 
well with observations. However, some dlff iclencies of this climate run 
are revealed: equatorial easterlies disappear, and zonal kinetic energy 

and momentum transport of stationary waves are too small in the model 
Northern Hemisphere (NH) . 

2 . The general features of wavenumber spectra of the model zonal and 
meridional kinetic energy for moving waves are consistent with observa- 
tion. The mode] zonal kinetic energy in wavenumber one of stationary 
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wave Is too small In the NH subtropics and equator. Furthermore, the 
zonal kirietlc energy of the model eastward moving waves in low wave- 
number regime is larger in NH and smaller in SH than observation. The 
meridional kinetic energy spectra of both hemispheres have smalLr 
values in low wavenumber regime and larger values In intermediate wave- 
number regime than observation. 

J. The frequency spectra of zonal and meridional kinetic energy of eastward 
moving waves show a band |5 'v 20 cycles (90 days) of large values in 
middle and high latitudes and peaks at 7 cycles (90 days) ^ 12.9 days. 

In addition, the frequency spectra have a slope of -2 "v -3 in a high 
frequency regime, rather than -1 as Kao finds. 

4. The model momentum transport in the wavenumber domain is mainly performed 
by the eastward moving and stationary waves of intermediate wavenumber in 
middle and high latitudes of the NH and only by the eastward moving waves of 
intermediate wavenumbers in middle and high latitudes of the SH. The soutli- 
ward momentum transport at the equator is due to the eastward moving 

waves of intermediate wavenumber. The momentum transport by the sta- 
tionary waves of low wavenumber is too weak in the model. In the fre- 
quency domain, the momentum transport is principaly carried by the east- 
ward moving waves of 5 20 cycles (90 days) ^ in both hemispheres, and 

1 'V 20 cycles (90 days) at 2*N. The momentum transport by the low- 
frequency waves is not significant in the model as found in observation. 

5. The model wavenumber-frequency spectra of zonal and meridional motion and 
momentum transport have a preferred band which is oriented from low wave- 
number and low frequency of the westward moving waves to high wavenumber 
and high frequency of eastward moving waves in middle and high latitudes. 


Thp preferred band of spectra tn the tropics extends from low wavenumber 
and low frequency to high wavenumber and low frequency. Furthermore, 
nuiximum spectral value of zonal motion appears In lower wavenumber than 
that of meridional motion. These features of model wavenumber^ f reqiiency 
spectra are consistent with the spectral studies of Kao, liayashi and Pratt 
Another Interesting feature of model spectra is that the maximum spectral 
value in NH occurs in lower wavenumber than in SH. 

The mean westerlies exist over the whole model globe at 200 mb In this 
climate experiment. The nonlinear interactions by the longitudinal 
convergence of kinetic energy flux (U1 and VI) and a combination of 
ageostrophic effect and dissipation (U4 + U5) and (VA + V5) are dominant 
and comparable. U1 and VI provides kinetic energy to the eastward 
moving waves, but extract kinetic energy from the westward moving waves. 
(llA -f U5) and (VA 4- V5) function oppositely as U1 and VI do. The primary 
interactions of kinetic energy in wavenumber-frequency domain are largely 
through the Interactions between mean zonal flow and the long and 
medium waves of low and medium frequencies. The significant role played 
by the stationary waves shown in Kao and Lee's does not appear In the 
model. 

The maintenance of momentum transport spectra by linear and non-linear 
processes in the NH is similar to that of kinetic energy spectra. However 
the functions of longitudinal convergence of momentum transport flux 
(VUl) and a combination of ageostrophic effect and dissipation (VUA + 

VU5) in SH is opposite to those in NH. At the equator, VUl supplies 
meridional flux of westerly momentum to both the eastward and westward 
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moving waves, but VUA VnJS extracts from both moving waves. This 
scheme is opposite to Kao and Lee's in direction. ’’’JiC primary inter- 
actions of meridional flux of westerly momentum is similar to those 
of kinetic energy. 
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Figure Captions 


Fig. 1 


Fig. 2 

Fig. 3 
Fig. 4 

Fig. 5 
Fig. 6 
Fig. 7 
Fig. 8 
Fig. 9 


. Latitudinal dlatrlbiition of zonal mean velocity, kinetic energy of 
zonal and meridional velocltiea of stationary and transient waves, 
meridional transport of westerly momentum at 200 mb. (a) GLAS model 
and (b) observation from various sources. The observations north of 
20*S come from KL, while south of 15*S are from E0LE study In which 
zonal mean velocity and kinetic energy are drawn from Morel and 
Desbols (1974), and momentum transport from Webster and Curtain 
(1974). The solid line denotes transient waves, while the dashed 
line represents stationary waves. 

. Wavenumber spectra of (a) zonal and (b) meridional kinetic energy 
of the GLAS model at 200 mb and various latitudes. 

. Same as Fig. 2, except for frequency spectra. 

. (a) Wavenumber co'^pectra of momentum transport associated with 

stationary and moving waves of the GLAS model at 200 mb and various 
latitudes, (b) Same as (a) except for frequency co-spectra. 

. Wavenumber-frequency spectra of zonal motion at 200 mb of the GLAS 
model. 

. Wavenumber- frequency spectra of meridional motion at 200 mb of the 
GLA.S model. 

. Wavenumber -frequency co-spectra of momentum transport at 200 mb of 
the GLAS model . 

. Schematic diagram for maintenance of spectral energy of moving waves 
at 200 mb of the GI.AS model. 

. Schematic diagram for maintenance of spectral momentum transport 
flux of moving waves at 200 mb of the GLAS model. 


Table 1. WavcnuMber Classification 


Wavenumber 

k 

Classification 

Symbol 

0 

zonal ly averaged flow 

0 

1-5 

long waves (small wavenumber) 

*‘s 

6-9 

medium waves 

*^m 

10-20 

short waves (large wavenumber) 



Table 2. Frequency Classification 


Frequency n 
(cycles per 90 days) 

Period 

(days) 

Classification 

(frequency) 

Direction of 
wave motion 

Symbol 

-90 to -31 

1-3 

high 

west to east 


-30 to -11 

3-9 

medium 

west to east 

^"m 

-10 to -1 

9-90 

low 

west to east 


0 


zero 

stationary 

0 

1 to 10 

9-90 

low 

east to west 

'"i 

11 to 30 

3-9 

medium 

east to west 

”m 

31 to 90 

1-3 

high 

east to west 

'"h 
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Lln««r And nonllnar contributions to snsrsy spoctrs In various oavsnuaboi-frs^usncy donsins at 
various lacitudss at 200 nb« 




ruu 

• U1 

♦ U2 

4 03 4 

(U44U5) 

BW 

• 

VI 

♦ 

V2 

a 

V3 a 

<V4aV5) 



15. as 

104.87 

♦ 15.50 

- 3.11 

. 

101.81 

28.55 


143.98 


4.28 


16.50 


94.03 


M*N 

45. 8S 

345.70 

- 81.87 

- 9.09 

- 

228.98 

8.80 


201.83 


4.18 

- 

8.60 

•• 

188.83 


1«*N 

18 . as 

40.99 

10.84 

4 0.20 

• 

33.38 

4.29 


15.88 

- 

0.20 

a 

3.13 ♦ 

14.00 


a** 

11.87 

5.79 

* 14.35 

4 0.01 

- 

8.47 

0.90 


1.07 

♦ 

0.43 

a 

U.Ol 

a 

0.61 

a 1 

li*S 

ao.38 

48.14 

4 90.72 

4 0.19 

- 

116.89 

2.91 


12.22 

• 

1.38 

a 

0.71 

- 

8.6/ 


38*S 

37.87 

243.63 

4 40.24 

4 1.72 
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Table’ 4. Linear and nonlinear contributions to northward momentum 
transport spectra in various wavenumber-frequency domains 
at various latitudes at 200 mb. 
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Abstract 


TVo experiments were performed previously by Shukla and Bangaru (1979^ 
to test the response of the general circulation model of the Goddard Laboratory 
for Atmospheric Sciences to a North Pacific sea surface temperature anomaly: 

(1) control run In which the climatological mean SST was used and (2) anomaly 
run In which an Invarient Pacific SST anomaly of 1976-1977 winter was imposed 
on the control run. It Is found that the ridges and troughs appearing In 
control run are amplified and become stationary in anomaly run. The purpose 
of this study is to examine the effect of SST anomaly on the long waves in 
terms of spectral analysis. 

The spectrally filtered HovmBller diagrams show that -avenumber 2 and I 
in the anomaly run become stationary around the climatological locations of 
observational study when the persistent blocking ridges were generated in this 
experiment. The constructive interference of wavenumber 2 and 3 in the anomaJy 
run forms two persistent blocking ridges: one in the west coast of North 

America and the other in western Europe. 

Spectral energetics analysis of the anomaly run during the appearance 
of persistent blocking ridges shows that A 2 and A-^, available potential energy 
of wavenumber 2 and 3, are supplied by the conversion from zonal available 
potential energy. K 2 and K^, kinetic energy of these two waves are maintained 
in different processes. K 2 is provided by C(A 2 » K 2 ) , conversion from A 2 to 
K 2 , while Kj is mainly furnished by -CCK^, Kr,), conversion from to Kj,, 
zonal kinetic energy. In other words, K 2 is supported by baroclinic energy 
conversion and is essentially supplied by barotropic energy conversion. 


The atmospheric circulation conditions in January 1963 and 1976-1977 
winter were similar to that in the anomaly run during the existence of per- 
sistent blocking ridge. The spectral energetics of these two periods of 
time are used to compare that of the anomaly run. It is concluded that the 
necessary condition to have blocking ridges formed is the constructive inter- 
ference of stationary long waves. The other condition for the formation of 
blocking ridges is how these long waves are amplified and maintained. The 
spectral energetics of blocking ridges shows that the amplification and 
maintenance of stationary long waves can be due to various physical processes. 
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1. Introduction 

The ocean haa long been known as an Important energy reservoir to the 
overlying atmospheric circulation which would be modified through tlie reii'ane 
of thermal energy (sensible and latent heat) from the ocean. Therefore, It 
Is Inferred that the appearance of the sea surface temperature (SST) anomaly 
may alter the atmospheric circulation pattern and generate anomalous weather. 
The observational studies in the past decades provided evidences of the large- 
scale feedback Interaction between the atmosphere and ocean. 

Bjerknes (1966, 1969) showed that the SST ano',aly over the tropical 
Pacific Ocean could cause the change of intensity of trade wind and Hadley 
circulation, and, in turn, might cause the anomalous atmospheric circulation 
In middle latitudes. Recently, Namlas (1978) indicated that the anomalous 
west-east SST gradients over the Pacific Ocean enhanced atmospheric baro- 
cllnlcity and also increased the upper level southerly flow. This anomalous 
southerly flow was not only to affect the local circulation pattern, but also 
the downstream circulation which might In turn affect the v;eather of Nortli 
America. Ratcliffe and Murray (1970) also found a feedback relationship 
between the Atlantic and western Europe. The appearance of negative SST 
anomaly over a wide area south of Newfoundland is associated with blocked 
atmospheric patterns the following month over northern and western Europe, 
while a positive SST anomaly in the same general area favor more progressive 
synoptic types. 

The observational studies described above indicates that SST anomalies 
may Influence vitally the climate dynamics and long-term forecast. To 
examine how the large-scale air-sea interaction affects the atmospheric circu- 
lation, one has to study the ocean-atmosphere system as an entity. Since the 
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whole atmosphere-ocean system is coo complicated, the general circulation 
model (GCM) has been used in the past decade as a quantitative tool to investi- 
gate the large-scale air-sea interaction. Moreover, Namias (1970) has noted 
that the SST anomalies, once formed, may persist for a long period of time, 
despite iterations with the overlying atmosphere. Therefore, the effects of 
the persistent SST anomalies on the atmospheric circulation are examined by 
inserting SST anomalies over the sea surface of a noninteractlve model. 

Rowntree (1972) used the Geophysical Fluid Dynamics Laboratory 9-layer 
hemispheric model to test the effect of the warm anomaly in the tropical east 
Pacific Ocean. The hypothesis proposed by Bjerknes (1966) was confirmed by 
Rowntree 's numerical experiment. The Goddard Institute for Space Study 
version of the 2-layer Mlntz-Arakawa model and multilevel version (Somerville 
ct al., 1974) were employed by Spar (1973a, b) and Spar and Atlas (1975) to 
test the effect of SST and SST anomalies on the long-term forecast. No signifi- 
cant impact has been found in their numerical experiments. 

The 5°-meslt version of the national Center for Atmcspherlc Research 
atmospheric model lias been used for a series of tests of the model response to 
SST anomalies. Houghton et al. (1974) tested the SST anomalies over the north 
At Lantic Ocean on the simulated climate. They found that the warm anomaly 
case showed good agreement with Ratcliff and Murray’s result. Chervln et 
al. (197b) imposed the dipolar SST anomalies on the Pacific. Their statistic 
analysis of 1.5 km temperature did not show downstream teleconnections 
between North Pacific SST anomalies and atmospheric temperature over the North 
American continent with high statistic significance. Kutzbach et al. (1977) 
made more extensive analysis for the experiments performed by Chervin et al. 
(1976) and also tested the effect of longitudinal location of SST anomalies. 


3 


In addition to the response of local circulation, wavellkt* perturb.it ionn uIho 
extended downstream from the region of SST anomalies. Chervin et al . (1980) 
made more experiments with various SST anomalies. Their experiment.*! showed 
that the midlatitude anomaly response involved near downstream wavelike 
changes and the su' iropic anomaly response was more confined to the vicinity 
of the SST anomaly. 

The observational study by Namias (1978) to examine the effect of the 
Pacific on the anomalous circulation pattern in 1976-77 winter indicated that 
the anomalous west-east SST gradients over the north Pacific was one of the 
major factors to cause a persistant open blocking ridge located in the west 
coast of North A.'nerlca. Following Namias' findings, two numerical experiments 
were performed to test the response of the Goddard Laboratory for Atmospheric 
Sciences (GLAS) atmospheric model to the SST anomalies over the north Pacific 
in the 1976-77 winter (Shukla and Bangaru, 1979). These two experiments are 
one using the climatological mean SST and the other imposed Invariant SST 
anomaly of 1976-77 winter. In addition to the climatological mean SST. The 
pronounced model response to the SST anomaly Is shown in the amplification 
of large-scale waves. The most remarkable aspect is the appearance of two 
persistent open blocking ridges located in tlie west coast of North America 
and western Europe. In order to understand how these two persistent blocking 
rldget are formed and maintained, our major effort in the present study is 
aimed at examining the alternation of the long-wave behavior due to SST 
anomalies and the spectral energetics of these long waves. 

A brief description of the GLAS model and numerical experiments is pro- 
vided In Section 2. Section 3 discusses the synoptic and spectral analysis 
of the experiments and the formation of persistent blocking ridges. Section 4 
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presents Che spectral energetics how the long waves which form the blocking 
ridges are maintained. A comparison of the model simulation to tho observa- 
tional analysis of the long waves and their spectral energetics Is made In 
Section 3. Section 6 offers the summary and discussion. 

2. Description of the GLAS general circulation model and experimental design 
Model 

The general circulation model used In this study Is the second-order 
version of the Goddard Laboratory for Atmospheric Sciences (GLAS) . This 
GCM was originally developed at Goddard Institute for Space Study (GISS) by 
Somerville et al. (1974) for a winter season simulation. A summer simulation 
with essentially the same GCM was carried out by Stone et al. (1977). The 
model was recently modified by Halem et al. (1978) at the GLAS in the model's 
longwave radiation parameterization, the surface boundary processes, and some 
numerical treatment of the time differential approximation. 

The model Is a global primitive-equation model which Is developed In 
the sigma vertical coordinate system (Phillips, 1957) with nine levels 
between the surface and the top at 10 mb. Horizontal grid spacing is 4® 
latitude and 5® longitude. A variable horizontal grid with coarser spacing 
near the pole Is introduced. In addition, the time-altering space-uncentered 
difference scheme of Arakawa (1972) was removed and a high order two-grid 
point filter developed by Shapiro (1970) was applied there. The time step is 
10 minutes in the model Integration. 

The model also includes a detailed parameterization of no explicit sub- 
grid scale processes. The condensation and cloud are generated by two different 
mechanisms: large-scale super saturation and parameterized subscale moist 
convection. Short-wave radiation Is calculated by a parameterization based 
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on detailed radiative calculariona, including multiple scattering. Long wave 
radiation fluxes are calculated using a new scheme developed by Wu (I960). 

These long- and short-wave radiation transfers were calculated at 5-liour 
intervals. The ground temperature are evaluated from the net heating and 
cooling at the surface by radiation, sensible heat and latent heat fluxes. 

The geographic distribution of the continent, topography, albedo, and ice and 
snow cover are specified. The soil moisture is predicted as treated in the 
UCLA model (CatcH, 1972). Sea surface temperatures adopt the monthly mean 
observed values provided by the NCAR (Washington and Thiel, 1970). These 
temperatures were updated dally by a linear Interpolation between the midpoints 
of two consecutive months. 

The verification of the CLAS atmospheric model simulation of winter and 
summer seasons can be found in Halem et al. (1978). 

Numerical Experiments 

In the winter of 1976-77, the severe weather occurred over North 
America, extreme cold dominated the eastern part of the United States and 
record drought affected the far west. In addition, abnormal warmth occurred 
over Alaska and much of the Canadian Arctic. Namias (1978) examined the 
possible factors to cause this severe winter weather. He found that the SST 
in the north Pacific was characterized by abnormally cold temperatures in its 
central and western parts with a warm pool off the west coast of the United 
States shown in Figure 1. In addition, the anomalous changes in SST are in 
phase vrlth the normal month- to-month changes. In other words, the normal 
SST gradient between 16t '-1A0“W was greatly strengthened in this winter. 

The enhanced west-east SST gradient over the Pacific increased barocllnicity 
in the overlying atmosphere and also increased the upper level southly oast 
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of 140®W. This flow helped s'eer storms northward to Alaska and produced a 
meridionally extensive ridge over the west coast of North America, and a 
trough in the east. Therefore, Namias suggested that the North Pacific SST 
anomalies might have been one of the multiple causes of the abnormally cold 
temperatures in eastern North America during the winter of 1976--77. 

A numerical experiment with the GLAS general circulation model was carried 
out to test Namias' hypothesis (Shukla and Bangaru, 1979). The experiments 
were performed by firstly integrating the GLAS model for 45 days with the 
climatological mean SST (Washington and Thiel, 1970) and the observed initial 
conditions for January 1, 1975. This integration is referred to as the 
c ontrol ru n. Ihe climatological SST field was then changed by adding a time 
invariant anomaly field and therefore, although the climatological SST varied 
with season, the imposed anomaly field remained constant with time. The 
distribution of the imposed SST anomalies is shown in Figure 1. Notice that 
the numerical values used in the experiments were exaggerated by considering 
tliera in °C rather than °F shown in Namias' study (1978). The model was inte- 
grated agaiii for 45 days and will be referred to as the anomaly run . Since 
this is not a coupled ocean-atmosphere model, SST anomalies were assumed to 
persist for the whole period of integration. This assumption is justified by 
the observed fact that SST anomalies do persist for a period of several months 
(Namias, 1970). The purpose of this numerical experiment is not to simulate 
the actual events of the 1976-77 winter but to examine the effects of similar 
SST anomalies on the model atmosphere. It should be noted that the atmospheric 


conditions in the control run correspond to the winter of 1975. 
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i. Synoptic and spectral analysis 
Synoptics 

The comparison of the synoptics between the control and anomaly runs will 
provide the first glance of the response of the model atmosphere to the SST 
anomalies shown in Figure 1. The examination of model dally charts (not shown) 
of both runs reveals that a two ridge-trough system always exist. Furtlier- 
more, the anomaly run generates two persistent blocking ridges in the west 
coast of North America and West Europe, in the first two weeks of model February, 
t.e, . the last two weeks of the anomaly run. In order to have an ove'/’vlew ot 
the synoptic difference between these two runs, the mean height fields for tlie 
first two weeks of model February at 300 and 830 mb are shown In Figure 2. 

Tlie two-week mean height fields display tlie circulation pattern of two 
ridge-trough system at 300 mb of both runs. The open troughs become close lows 
at 850 mb. Two ridges exist along the west coast of North America and western 
Europe at both 300 and 830 mb. The ridges of the anomaly run are amplified 
very remarkably compared to the control run, especially at 500 mb. Tlie 500 mb 
trough over the Atlantic is deepened very much In the anomaly run. The 300 mb 
trough over the Pacific is very broad. However, ’•.he 850 mb close low, which 
is 500 mb open trough, over the Pacific also demonstrates its deepening in the 
anomaly run. It is also interesting to note that both the trough and ridge 
lines Kit from southeast to northwest. In addition, the ridge of tlie anomaly 
run over western Europe Is shifted westward at 500 mb and eastward at 500 mb 
with respect to that of the control run. This contrast between both runs is 
especially pronounced north of bO'N. 
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TIu’ ikH'pcnlnR of the low (or trouglO over the I’acifle and the Intenslfl- 
cat Ion of the ridge over the west coast of North Air.erica can be understood due 
to tl>e enhancement of the model barocHnlclty by tlie imposed west-east SST 
anomalv gradient and also possibly by the Interaction between the positive SST 
anomalies over the east North Pacific and the Rockys. !o other worde, the 
persistent b locking ridge appearing on the west coast of North America can be 
intuitively uviderstooil as the near dowistream teleconnection. Furthermore, 
t tu> amplification of the ridge over western F.urope and the deepening of the 
t rnugli (or low) over tlte Atlantic naturally demonstrates ttuU the (IbAR model 
atmospliere also respond to tl\e Pacific SST anomaly wltl> a far downstream 
t I'leconnect Ion . 

Spectral Analysis 

lU)vmoller's diagram, i.e. longltude-t line trough-ridge system, (llovmoller, 
ld/t‘1) has li'ug been usi'd to examine tlie amplification and in-opagat ion of 
at mosnlu'C ic long waves. Recentlv, Somerville (1^)8(V) suggested tliai the 
spectrallv filti’red (SF) IU>vmol ler diagram is a more effective mean to isolate 
t lu' ii\teresling leatures of long waves and to examine tlie four-dlmenslonal 
development of long waves. 

Till' comparison ol the height fields between both runs reveals clearly 
that the most distinguishable difference is In the large-scale features. In 
other words, the major response of the GLAS model atmosphere to the SST 
anomalies over the North Pacific show.s in the long waves. For this reason, 
it would be instructive to use the SF llovmoller diagram to illustrate the 
spectral response of the model atmosphere to tb.e SST anomalies. Furthermore, 
')0“N, which Is located in the middle of the north-south extent of the ridges. 
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is chosen to depict the SF Huvraitfller diagrams. Figure 3 shows these diagrams 
of three different wave groups: wavenumber 1-15, wavenumber 1-4, and wave- 

number 5-15, at 500 mb. The first wave group is regarded as the total wave 
disturbances, the second wave group as the long waves, and the third wave 
group as short waves. 

The SF Hovmdller diagrams of wavenumber 1-15 for both the anomaly and 
control runs. Figures 3a and d, show that two ridge-trough systems exist over 
most of the simulation period. The trough-ridge system in the anomaly run 

stationar y , while that Is eastward propagating in the conf :ol run. Comparison 
between Figures 3a and d also reveals that the trough-ridge systems in both 
runs in the third week of January are still similar to some extent. Tills 
similarity of the SF Hovmdller diagrams for total wave disturbances may indi- 
cate that the GLAS model has memory of the initial height field and the model 
response to the Pacific SST anomalies does not become significant in the first 
three weeks of the model simulation. 

In order to examine the role played by the long and short wave regimes, 
the SF Hovmdller diagrams of wavenumber 1-4 and 5-15 are also displayed in 
Figures 3b, c, e and f, respectively. The comparison between the long-wave 
SF Hovmoller diagrams and the total-wave SF Hovmdller diagrams shows very 
clearly that the large-scale pattern of the ridge-trough system in both run.s 
are mostly explained by the long-wave regime. The contribution from the short- 
wave regime is the fine structure of the ridge-trough system. The long-wave 
SF Hovmoller diagrams show that the persistent blocking ridges in the anomaly 
run are mainly constructed by the stationary long waves. 

The short-wave SF Hovmdller diagrams characterize the eastward propagation 
of the short waves over most of the longitudes and time. However, both the 
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control and anomaly runs have stationary ridges located at 30®E in the first 
two weeks ut February* In addition, the troughs and ridges of the short- 
wave regime in the anomaly run between 150®E and 120°W are quasi-stationary . 
Recall that the SST anomalies are imposed in the North Pacific with the 
boundary of positive and negative anomalies around 140°W. It is understood 
tli.it the strong west-east SST gradient would enhance the baroclinlcity of the 
overlying atmospliere. This may, in turn, increase the cyclogenesis in the 
enhanced baroclinlc zone. In any event, the quasi-stat ionar ity of the short- 
wave ridges and troughs may indicate the statlonarity of forcing from the 
Pacific SST anomalies. This argument is consistent witli previous experiments 
using the NEAR atmospheric model by Chervin et al. (1976), Kutzbacli et al. 
(1977) and Chervin et al . (1980). The Intensification of the ridge located 
at 30“E is a very surprising downstream teleconnection in the G1.AS atmospheric 
model . 

To further our understanding of tlie long-wave response to the Pacific SST 
anomalies, the SF Hovmollcr diagriuns for each Individual Jong wave from wave- 
luimber 1 to 4 are dejiii-ted in Figure 3. The most pronounced difference 
between two runs are wavenumbers 2 and 3. These two waves "lock In" the 
anom.ily run. In other words, the west-east Pacific SST anomalies force these 
two waves to be stationary. It is also interesting to note that the super- 
imposition of wavenumbers 2 and 3, Figure 3e , provides a reasonable expla- 
nation of the double persistent blocking ridges, shown in Figures 2a or b. 

Based upon the occurrence of bjecking observed by Rex (1958) and the 
latitudinal distribution of planetary-wave phase at 500 mb analyzed by Eliasen 
(1958), Austin (1980) claims that the formation of blocking can be classified 
in the following ways due to the constructive interference and amplication of 


planetary waves: 
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Wavenumber 1 Wavenumber 2 


Wavenumber 3 


Atlantic Blocking large 
Pacific Blocking small 
Double Blocking large 


large 

large 

large 


small 

large 

large 


Although the anomaly run does not generate the Rex-type blocking, the 


amplification of two persistent blocking ridges In the first two weeks of 


February seems to be similar to double blocking; one ridge Is located between 
the eastern North Pacific off the west coast of the North American continent, 


and the other is situated between the eastern North Atlantic and western 


Europe . 


It would be interesting to compare the spectral analysis of the anomaly 
run with Austin's compilation of observational studies. Wavenumber 1 in the 
anomaly run Is shifted about 50“ westward with respect to Eliasen's result. 
Wavenumber 2 of the anomaly run has about the same location as observation, 
witile, wavenumber ! of the anomaly run is shifted eastward about 10“ . The 
SF Hovmbller diagram of wavenumber (2+3), Figure 4e , shows clearly that the 
eastward shift of wavenumber 3 makes the formation of double blocking ridges 
only by wavenumbers 2 and 3 possible. The amplification and westward shift of 
wavenumber 1 helps the filling of the trough between two persistent blocking 
ridges formed by wavenumbers 2 and 3, and the intensification of the Atlantic 
blocking ridge. 

Tung and Llndzen (1979) proposed resonance theory to explain the develop- 
ment of atmospheric blocking. "To be resonant, the wind condition in the 
atmosphere has to be such that the phase speed of the free traveling wave is 
reduced to zero, so that the wave becomes stationary with respect to the surface 
of the Earth and hence, also to the topographic forcing and land-sea differential 
heating." This is Tung and Llndzen's "stationary condition". The discussion 
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above clearly shows chat the insertion of the Pacific SST anomalies make 
wavenumbers 2 and 3 "lock in". This lock-in effect enable wavenumbers 2 and 
3 In tlie anomaly run satisfy the necessary condition for the occurrance of 
Clie persistent blocking ridges. Furthermore, one has to face another problem, 
chat is the amplification of the stationary waves. Austin (1980) speculated 
that the imposed SST anomalies enhance the theniial forcing which, in turn, 
amplifies stationary waves to form blocking. In order to understand this 
problem, we shall apply the spectral energetics to examine how long waves are 
amplified after becoming stationary in Section 4. 

It has been sliown that wavenumbers 2 and 3 in the anomaly run are so vital 
to determine the formation of persistent blocking ridges. It would be of 
Intcri'st and informative tt) examine t lie synoptic structure of tliese two waves, 
iiie persistent blocking lidges obtain their maximum intensity on February 4. 

It would be illustrative to show the spectrally filtered synoptic charts of 
wavenumbers 2, 3 and combination of them at 12 GMT 4 Feb for height and 
temperature at 500 mb, Figure 3. 

The highs of wavenumber 2 located at Alaska Bav and north Central Europe, 
while those ot wavenumber 3 are anchored at the west coast of Canada, the 
North Sea and Siberia. The combination of these two waves forms two highs 
at the west coast of Canada and the North Sea with a deep low located at 
i,abrador and another broad low extending from north Africa to the northwestern 
Pacific. The combination of temperature wavenumbers 2 and 3 is very similar 
to the height except the waniith centered at the Alaska Bay is shifted westward 
with respect to the high tu height field. The comparison of individual waves 
of height and temperature shows that wavenumber 2 has a significant phase lag 
between height and temperature and wavenumber 3 does not. In fact, the 
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observations of daily spectral analysis reveal that wavenumber 2 height le.uls 
wavenumber 2 temperature from 10* to 30® longitude. On the contrary, heli'.ht 
and temperature of wavenumber 3 does not show very significant phase lag on 
the basis of dally variation. This indicates that wavenumber 2 is barc'cllnlc 
and wavenumber 3 is relatively barotroplc. The remarkable difference of 
syno|)tlc structure between these two waves must be related to the difference of 
physical processes to develop and maintain them. 

4. Energetics 
E nergetics Scheme 

In order to examine the development and maintenance of the blocking ridges 
in terms of spectral analysis, Saltzman's (1957) spectral energetics scheme Is 
employed to a latitudinal zonal between 26®N and 86*N. Following Saltzman's 
formulation and Steinberg et al.'s (1971) expression, spectral energy equations 
In symbolic form can be written as: 

N 

I C(K^, K^) + C(A2, K^) + F(K2> " IKK^) , (1) 

n=l 

- C(K^, K^) + CK(n|m,0 + C(A^, K^) + F(K^) - D(K^) , (2) 

N 

- I ClA^, A^) - C(A2, Kj,) -t F(Aj,) + fi(A^) , (3) 

n*l 

C(A , k) - C(A , K ) C(nlm,«) + F(A„) + G(A„) , 

^ n n IX n n 




dt 


dK 


dt 


dAj, 

■dF 


dA 
n 

dt 


(4) 




when* K and A represent kinetic and available potential energy. The subscripts 

Z and n denote zonal and wave components. The rate of change of various energy 

components on the left-hand side is contributed by the various energetic 

processes on the right-hand side. C(P,Q) represents the energy conversion 

N N 

f£ *m P to Q. Note that ** ^ ^n^ " • 

CX(n/ni,«) is the gain rate to wavenumber n of quantity x due to triplet inter- 
actions between all possible combinations of m and it with n. F(x) gives tlie 
contribution due to the flux across the lateral boundary of the open region. 
C(A) and D(K) are generation of available potential energy and dissipation of 
klnctie energy, respectively. 


N N 

Equations for eddy energies, = 'I K and A^, ■= ^ A^, can be obtained 

n*l ' n*l 

by summing equations (2) and (4) from n = 1 to N, 


dK 

== - C(Kj^, Kj,) C(Ap, Kg) + F(Kp) - D(K^), 


(5) 


dA 


( 6 ) 


The traditional Lorenz energy cycle consists of equations (1), (3), (5) and 
(6). A schematic energy diagram is shown in the appendix to Illustrate the 
maintenance of wavenumber n energetics. 
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Overview of model energetlcB during blocking sltuatlonB 

Time variations of various cmergles, A^, and Kg, and various 
energy conversions, C(Ag, Ag), C(Ag, Kg) and C(Kg, Kg) from January 16 to 
February lA In the anomaly run are depicted In Figure 5. Lei; us recall the 
amplification of blocking ridges located In the west coast of North America 
and western Europe in the anomaly n » covers the first two weeks of February. 
Observation of Figure 6 reveals that A,, and K„ start to decline on January 30, 
wiiile Ag and Kg increase. Then, A^. id Kg, especially the latter, decrease 
around February 7, and both increase about February 10 again. Notice that 
both Ag and Kg increase. In general, when Ag and Kg decrease, or vice versa. 

The energy conversions, C(Ag, Ag) and C(Ag, Kg) are always positive. 

In addition, their time variations go hand-in-hand and are also consistent with 
the time variations of A„ and K_. ft Is surprising to note that C(K„, K„) 

Is negative over most of the time period In the anomaly run. This indicates 
that the eddy disturbances are barotropically maintained by the energy conver- 
sion from the zonal flow. 

Lejenas (1977) compiled the time variations of A_, A„, K„ and K_ from 

L 't L b 

the energetics study of atmospheric blocking by Winston and Kreuger (1961) and 
Miyakoda (1963) and energy conversions by Paulin (1968). Lejenas sketched a 
characteristic picture of the energy changes and energy conversions during a 
blocking situation in Figure 7. We may not be able to expect that the energetics 
for the development of the persistent blocking ridges in the anomaly run is 
the same as the atmospheric cases analyzed by either Winston and Kreuger or 
Miyakoda. However, the energetics during the development of model blocking 
ridges depicted in Figure 4 shows a very striking similarity to LejenSs’ 
schematic picture. 


I)jMiamlcal t f an spo rt s of wavenumbers 2 and 3 during blocking 

It was observed In Section 3 that the life period of the blocking ridges 
generated In the anomaly run was about two weeks. In addition to examine hew the 
blocking ridges are formed, we have to face how the blocking ridges are main- 
tained for such a long period of time. The spectral analysis in Section 3 
shows chat the formation of the blocking ridges are essentially due to the 
cnnsl ructlve Interference of wavenumbers 2 and 3. On the other hand, the 
spectral energetics of these two waves would shed light on the maintenance of 
blocking ridges. 

Before we Investigate the spectral energetics of wavenumbers 2 and 3 in 
detail, the examination of the dynamic transport activities would provide some 
Informative property of these two waves. Figure 8 depicts the time variation 
of xonal flow and temperature at 500 mb for the anomaly ruuo The zonal flow 
between 2b°N-62“N starts to decrease when the blocking ridges develop at the 
end of January. The zonal flow Increases again after about A days. An 
Interesting feature about the zonal flow Is the appearance cf easterlies in high 
latitudes (60“ - 70°N) from February 3-9. This Is obviously caused by the warm- 
ing at high latitudes shown in Figure 8b. The warming at high latitudes 
indicates that the blocking ridges channel the warm air to high latitudes. 

The time variations of momentum and sensible heat transports* for waves 
2 and 3, uv(2), uv(3), vT(2) and v^(3), are displayed in Figures 9 and 10, 
respectively. It can be seen clearly that the transport activities by wave- 
numbers 2 and 3 mainly occur in high latitudes. This is consistent with the 

*vf(n) and uv(n) will be used hereafter to denote the sensible heat and 
momentum transport by wavenumber n. The combining contry>ution from two waves, 
e.g. 2 and 3, will be denoted by vT(2+3) and uv(2+3). ( ) = zonal average. 
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obBervational specectral analysis of these two transports by WMn-Neilsen 
et al. (1964). It Is interesting to note that uv(2) aitd uv(3) (vT(2) and vT('J)) 
are, in general, southward (northward) in high latitudes. Surprisingly, uv(3) 
becomes northward from midlatitudes to high latitudes during the mature stage 
of blocking ridges. vT(2) is very significant northward transport during the 
existence of blocking rl.lges, but vf(3) is significantly southward over sonit- 
period during the developing stage of blocking ridges. This may manifest that 
vT(2) Is Important fac'or to cause the warming in high latitude shown in 
Figure 8. 

The blocking ridges shown in Figure 2 reach their maximum Intensity on 
February 4, The geographic distribution of uv, vT and u)T of wavenumbers 2 and 
3, and sum of these two waves at 500 mb on February 4 are also shown in Figure 
11. Observation of these figures reveals that uv(3) is more significant than 
uv(2), while vT(2) and u>T(2) are larger than vT(3) and o)T(3). This contrast of 
transport properties between wavenumbers 2 and 3 Is consistent with the synoptic 
analysis that waveniunber 2 is baroclinic and wavenumber 3 Is relatively barotropic. 
uv(2+3), vT(2+3) and u)T(2+3) are significant In the western hemisphere, it is 
interesting to note that vT(2+3) and oaT(2+3) show three areas of large value: 
eastern Pacific, North America, and eastern Atlantic. The superimposition of 
height and temperature geographic distributions for the sum of wavenuiabers 
2 and 3, shown in Figure 5, explain well the distribution of vT(2+3). Positive 
vT(2+3) at eastern Pacific and eastern Atlantic is due to northward warm air 
advectlon, while positive vf(2+3) at North America is because the southward 
cold air advection. -u)T represents the vertical transport of sensible heat 
which usually is in phase with vT. In addition, -o)T also denotes the conversion 
between Aj. and K^., namely, the release of available potential energy. The 
"lock-in" effect by imposing Pacific SST anomalies in the anomaly run seems 
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Lu increase the release of to enhance stationary wavenumbers 2 and 3 and, 

In turn, to amplitv the intensity of blocking ridges. 

Kutzbach et al. (1977) imposed a dipole SST anomaly on the Pacific of 

the NCAR model with warm part to the east and cold part to the west. The 

spacing between the centers of warm and cold SST anomaly is about 45° of 
longitude wlUch is about wavelength of wavenumber 4. The model temperature 
clianges have tlie greatest magnitude and vertical extent over the warm ocean 
anomaly, rather than a wavenumber 4 response. Kutzbach et al., therefore, 
speculated that the warm anomaly and Its longitudinal position with respect 
lo the continents may be of primary Importance to describe the pattern of 
mid-troposphere response. The geographic distribution of vT(24-3) and mT(24-3) 
in the present study may support their conjecture. 

Spectral enerjgetlcs of wavemimbej .2^_ajfd_3 

A and K of various long waves In the anomaly run as a function of time 

are sliown in I’lg. 9. The inclusion of wavenumber 1 and 4 is to provide a 

I'omparfson i)f energy contents in different long waves. Although Aj Is dominant 
over some period In f'ebruary, the SF Hovmdller diagrams show tiiut the formation 
of the persistent blocking ridges is mainly due to the constructive inter- 
ference of WlIvi number 2 and 3. Note that both A and K of wavenumber 2 and 

n n 

3 Increase from the end of January when the persistent blocking ridges shown 
in Fig. 2 start to develop. Furtbeniiore , our spectral energetics analysis 
Indicates that the major contributions to various energy conversions of Lorenz 
energy cycle in the. long wave regimes are from wavenumber 2 and 3 in the anomaly 
run. 

Let us first deal with the energetics of wavenumber 2 with the aid of 
spectral energy diagram displayed in the appendix. A,, is maintained by 
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C(A^, A2>, C(A2» K2> and CA( 2 jin,)l). Both C(A^, A2) and C(A2, K,) are 
poalLlva over the time period concerned In our analysis. Furthermore, the 
time variation of C(Ag, A2) and C(A2, K2) are In concert with that of A2. 

This In-phase variation Indicates that A2 Is provided by C(A^, A2)« but Is 
consumed by C(A2* K,,) . The nonlinear Interaction CA( 2 |m ,0 actually also 
extracts energy out of A2 over most of tne blocking period. 

K2 Is virtually maintained by C(A2» K2), C(K2» K^,) and CK( 2 jm,it). 

C(A.,, K2) was found to be always positive. C(K2* K^) is positive over most of 
the time except in the second week of February. CK( 2 |m ,0 is generally 
negative except when K2) drops significantly during February 8 - 9 . It 

is obvious that C(A2, K2) maintains K2» l^(K 2 » ^ 7 } CK( 2 |m,l) diminish 

K2. 

C(A2» A2) and C(A.^, K2) Involve the barodinlc processes of horizontal 
and vertical sensible heat transport, respectively, while 0(K2» K^) is due 
to the barotroplc process of momentum transport. The maintenance of A2 and 
K2 are provided by C(A2, A2) and C(A2» ’ *^®spectlvely . Recall that the 

synoptic analysis in Section 3 showed that wavf’number 2 temperature alwavs 
has a phase lag behind wavenumber 2 height. These observational evidences 
indicate that the development and maintenance of wavenumber 2 is duo to the 
bar o clinlc process . 

The time variations of various energy conversions to maintain A^ and 
are depicted in Figure 11 . It was shown in Figure 9 that A^ does not increase 
very much compared with A^ and (<2 during the development of blocking ridges, 
but increases very significantly. Let us follow the argument for the 
maintenance of A2 and K2 to examine the energetics of wavenumber 3 . C(A2, A^) 
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and CA(3|m, il) are major channels to support However, the small value of 

C(A^, A^) compared with explain the reason why A^ is not increased 

very much during the development of the blocking ridges. Furthermore, the magni- 
tude of CCAj, K^) is not large either. It may be because CCA^, A^) is small. 

is essentially supplied by -C(K^, , the barotropic conversion releasing 

the available zonal kinetic energy, rather than C(A 2 , K 3 ) , the baroclinic con- 
version to release the eddy available potential energy. It is very clear that 
the physical processes tc maintain the energetics of is different from those 
to K 2 . As we discussed previously in Section 3 that wavenumber-3 height and 
temperature do not possess any significant phase lag between them. In addition, 
the momentum transport of wavenumber 3 becomes very significant during the 
development of blocking ridges. In any event, these arguments provide evidences 
that wavenumber 3 is mainly developed and maintained baro tropically during the 
existence of blocking tidges in the anomaly run. 

The contrast of physical processes to develop and maintain wavenumber 2 
and 3 during the existence of persistent blocking ridges in the anomaly run 
is very apparent. In fact, the energetics of the blocking ridges occurred in 
January 1963 are similar to that generated in the anomaly run to a great extent. 
The spectral energetics study for this peculiar month, which will be discussed in 
further detail in the next section, by Wlin-Nielsen et al. (1964), Wiin-Nielsen 
(1964) and Murakami and Tomatsu (L965) showed that wavenumber 2 and 3 were 
dominant in this month. Wavenumber 2 were maintained mainly by baroclinic pro- 
cesses, while wavenumber 3 were supported by barotropic processes. 

5. Comparison with observational studies 

The model atmosphere is not a real one and the model response to the 
SST anomalies may not be the same as the response of the real atmosphere. 
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Nevertheless, the comparlwon u£ the present analysis of the anomaly run 
with some observational studies may shed some light on the physical mechanism 
which causes the persistent blocking ridges generated in the anomaly run. 
Under this guidance, we select January 1963 and January-February 1977. 

The reasons to select these two winters are as follows. 

It was pointed out by Tung and Lindzen (1979) that the United States 
experienced very similar abnormal circulation patterns and severe winter 
conditions In January 1963 and January 1977; a persistent and enormously 
amplified ridge existed in the eastern Pacific. Furthermore, the spectral 
energetics of the anomaly run during the existence of tl>e persistent blocking 
ridges is similar to that of January 1963 analyzed by Wlln-Nlelsen et al. 
(1964) and Murakami and Tomatsu (1967) to a great extent. Examining the 
large-scale interaction between atmosphere and North Pacific, Namlas (1978) 
suggested that the North Pacific SST anomalies in 1976-1977 winter could be 
the major factor to maintain the quasi-statlonary amplified long waves over 
the North Pacific and North America. The anomaly run in the current study 
adopts the North Pacific SST anomalies compiled by Namias (1978). Of courue, 
we do not expect these SST anomalies can force the modei to generate the 
same circulation condition as the real atmosphere in 1976-1977 winter. 
However, it would be Interesting to see whether the model response to tlie 
Pacific SST anomalies in the same direction as Namlas' suggestion of the SST 
anomalies on the atmospheric circulation. 

(a) January, 1963 
Synoptics 

It was observed that the existence of an unusual type of January 1963 
atmospheric circulation which had an extreme amplification (O'Connor, 1963). 


riu> aviMMgo ivlrculat Ion of 700 mb Is shown in Figure 14. The pronounced 
liMluroH of thv at moH pile tic circulation in this month are enormouHly 
amplified ridge in the eaBlern Faclfle, a blocking higba of Rex type pre- 
vailed through the month between Iceland and IJrltain, ai.d mother strong 
ridge over llie central iVsia. These three ridges were accompanied with the 
abnormally amplified eastern North American trough, deep A/.ores low, 

Kiiropean trough and tlie pacific low extending from the west Pacific to 
Hawaii. Tlie cold artlc air was brought to the nortlu ast North America, 

Kurope and Far Fast, .uul caused very severe cold weather over tliese three 
aiv-as. The unusual warmth was occurred over Alaska and Creeniand. In 
addition, most ot the trough and ridge lines tilted from south-east to north- 
west, This I’irculatlon condition Is similar to the existence of the per- 
sistent blocking rldg.es in the anomaly run. 

i'hierge tics 

The I'nergetics analysis of .lanuary 19b 1 discussed in the following based 
upon studies ot Wi in-Nielsim (1964'), Wlin-Nlel' cn al. (1:U>4) and Murakami 
and Tomatsu (1967). The horenz energy cycle of this monili is slunm in 

_2 

Figure lac. Tlu' directions ot energy conversions, t’(A^, A,,) = 6.4 1 Watt m 
and C(Aj,, Kj,) - 2.96 Watt m are the same as convi'iit tonal studies (Lorenz, 
1967). Surpr i s i ngl V , C(K|,.. K^) ot tills month has a large nej>ativi' value 
i)f -1.29 Watt m “■ which indicates thiit the eddies were barot rop i ca 1 1 y unstable 
in this month. Tlu> time variations of various energy component.s in danuarv 
196 3 are depicted in Figure 15a and b. Note that A^, and K^, Increase (decrease) 
when A,, and K„ decre.ise. The overall picture of time variation of January 
1963 atmospheric energetics shown in Figure 15 is similar to Lejeniis' during 
a blocking situation described in Figure 7. 
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The spectral energetics of January 1963 is displayed In Figure 16. 

The spectra of A^, Figure 16a, and Figure 16b, show that wavenumber 2 

and 3 have the largest values, although is comparable to A 2 or A^. This 

may Indicate that the atmospheric circulation of this month is mainly 

described by wavenumber 2 and 3. It may also manifest that the persistent 

blocking ridges in January were due to the constructive interference of these 

waves. C(Aj^, K^) is extracted from the result of Murakami and Tomatsu (1965) 

which was integrated over 30“N-70“N and 1000-500 mb. In any rate, it still gives 

some Indication of over the whole layer of atmosphere. C(Aj,» » 

Figure 16c, and C(A^^, , Figure 16e , of wavenumber 2 and 3 are positive and 

dominant, while CCK^, , shown in Figure 16d, has very large negative value 

-2 

ol -1.14 Watt m . Comparing Figure I6d and Figure 15c, we can see that 
-eCK^, K^) explains most part of -C(Kp, K^). This shows that wavenumber 3 
Is barotropical ly unstable. These discussions above Indicate that A,, and A^ 
are supported bv C(A^, A 2 ) and C(A^, A,^), respectively . Note that A,,) 

> 0(Ar,, A^) may explain why Ay > A^. Furthermore, K 2 is maintained by 
0{Ay, K 2 ), baroclluic energy conversion, while is sponsored by both C(A^, K.^) 
and -C(K^, K 2 )* other words, is supplied by a very significant baro- 

tropic conversion in addition to baroclinic conversion. The energetics of 
wavenumber 2 and 3 in January 1963 is very similar to that in the anomaly run. 

Lot us recall that wavenumber 2 and 3 "lock in" by imposing the Pacific 
SST anomalies of 1976-1977 winter in the anomaly run. The spectral energetics 
of wavenumber 2 and 3 in January 1963 and the anomaly run is simx!.';;,% We may 


speculate that tlie mechanisms to cause the "lock-in" of these two waves in 
these two cases may be tlie same. However, the pattern of Pacific SST anomaly 
in the 1962-1963 winter was an extensive and persistent warmth over central 
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and eastern North Pacific (Namias, 1963) which Is not the same as that of 
1976-1977 winter. 


( B ) January - February 1977 
Synop ti cs and spectral analysis 

The main features of atmospheric circulation in 1976-77 winter were very 
persistent through the wliole season (Taubensee, 1977; Wagner, 1977; Dickson, 
1977). The mean 700 mb circulation of this winter is shown in Figure 7 
(Namias, 1978). Deep troughs persisted over the central Pacific ocean and 
eastern North American while a strong persistent open blocking ridge located 
near the west coast of North America. Therefore, the warm air was channeled 
to Alaska by the southerly to the west of the ridge and the cold artic air 
accompanying with heavy snows was brought to the northeast and droughts also 
occurred in the west, in addition, some blockings occurred over the North 
Atlantic and the pole in January and February. 

Tlie spectral analysis of wave motions are illustrated by the 500 mb SF 
Hovmoller diagrams at 50°N, Figure 18. The total-wave SF Hovmbller diagram 
presents two very persistent blocking ridges: a very strcng narrow blocking 

ridge between 90‘’W-150°W and a broad weak blocking ridge between 30“W-90°E. 

The long-wave SF HovmOller diagram explain the major feature of the blocking 
ridges as we have seen in the anomaly run. The strong narrow blocking ridge 
between 90°W-150°W existed since the first week of January 1977 and decayed 
at the end of the third week of February. Thb short-wave SF Hovmoller diagram 
show that the small waves have strong eastward propagation. However, exception 
occurs between 120®W and 120°E, i.e. over the Pacific, where the ridges and 
troughs composed of short waves becomes relatively stationary. Namias (1978) 


claimed that the significant west-east Pacific SST gradient enhanced the 
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cyclogensis and the upper level southerly which steered the storm northward. 

The quasl-statlonarlty of short-wave ridges and troughs may manifest that 
the cyclogensis occurs In the same areas very often. This feature Is also 
found In the anomaly run. 

The SF HovmfJller diagrams of Individual waves are also displayed In 
Figure 15d-h. All long waves (n»l-4) , during January-February 1977 at 
50“N were quasi-statlonary . The two-month mean locations of these long- 
wave ridges are very close to Ellasen’s (1958) values which are displayed in 
Figure 1 of Austin (1980). It is Interesting to note that the very strong 
narrow blocking ridge decayed when wavenumber 3 and 4 became moving and 

I 

j 

k weak. 

\ 

I Let us compare the SF Hovmoller diagrams of individual long waves in the 

anomaly run. Figure 3, and those during January-February 1977. The ridge of 
model wavenumber 1 is located around bO^W, when this wave becomes stationary 
in the anomaly run. The ridge of observational wavenumber 1 is located 

I 

around 0® of longitude. This discrepancy may be due to the difference of 
seasonal forcing between the model simulation and observation. When wave- 
number 2 and 3 become stationary in the anomaly run the ridge positions of 
these two waves are very close to Eliasen's observation. Wavenumber 4 is 
essentially eastward moving wave in the anomaly run. The effect of the Pacific 
j SST anomalies make wavenumber 2 and 3 stationary in the positions as these 

two waves found either in January-February 1977 or in Eliasen’s (1958) 

* 5 . 

4 studv. 

I 

1 
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Kne t ics * 

Our c'uergetlos computations for January-February 1977 reveal that time 
varlallouH of Aj,, and were not very significant, and A^, and Kj, increased some- 
what in February (not shown). The slight enhancement ol A., and K^, Is due to the 
intensity deireasi- lO the persistent trough and ridge in this winter (Dickson, 
1977). Figure 18 shows that an Atlantic blocking existed in the first ten 
days of January. Time variations of various energies during this period 
of time are the same as Lejeniis' scliematic diagram. Figure 5. 

I.orenit energy cycle of these two months is depicted in Figure 19. The 
directions of various energy conversions are the same as conventional studies. 

A remarkable difference between the energy cycle of January-February 1977 and 
that t)l the anomaly run is C(Kj^,, Kj-) ' negative in the anomaly 

run, while it is positive in January-l’ebruary 1977. In other words, the 
eddies of this winter is barotropical ly stable, Namias (1978) suggests that 
the seasonal forcing i' another factor to decide the atmospheric circulation 
pattern. The contrast of C(Kp, K^,) between January-February 1977 and the 
anomaly run seems to indicate again that the seasonal forcing of the GLAS 
model is difterent from 1976-1977 winter, in addition, Rosen and Salstein 
(1980) pointed out that 1976-1977 winter has strong C(Kp, R^) which is due 
to the highly amplified nature of atmospheric waves. 

The spectral energetics of January-February 1977 is displayed In Figure 
20. A^ is somewhat smaller than Aj^ and A 2 . but is comparable to and 

*The computations use the observational data analyzed by the National 
Meteorological Center (NMC) and archived at the National Center for 
Atmospheric Research (NCAR) . The region of integration covers 30°- 
80°N. The u> fields,, which were used to calculate C(A , Rp) * were 
evaluated by the scheme proposed by Chen et al. (1981J. 
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blocking ridge In the west coast of North America is mainly due to the con- 
structive interference of wavenumber 2 and 3. In order to compare with tlie 
energetics of anomaly run, we shall focus on these two waves. C(A^, A 2 ) and 
CCA^, Aj) dominate C(A^, Aj,) and are the sources to supply A 2 and A^, 
respectively. K 2 and are maintained by C(A 2 , ^ 2 ) and C(A^, K^) , respec- 
tively. The most pronounced difference of spectral energetics between these 
two waves is that has a large energy sink, C(K^, K^) which explains 45% 
of C(Kg, K^) . Note that C(K 2 , K^) is small compared to C(K,j, K^) • The 
comparison of C(K^, K^) between this winter and the anomaly run shows that 
wavenumber 3 is barotropicaily stable in January-February 1977, and barotroni- 
cally unstable in the anomaly run. In other words, wavenumber 3 has opposite 
physical functions in these two atmospheric circulations. 

6. Concluding remarks an d dis cussion 

The general circulation model of the Goddard Laboratory for Atmospheric 
Sciences was used to test tlie model response to a Pacific sea surface 
temperature (SST) anomaly. Two experiments were performed: (1) control run 

in which the climatological mean SST compiled by Washington and Thiel (1971) 
was used, and (2) anomaly run in which an invarient Pacific SST anomaly of 
1976-1977 winter compiled by Namias (1978) was Imposed on the control run. 

The ridges and troughs appearing in control run are amplified and become 
stationary in anomaly run. The most pronounced difference between these two 
runs is that the Pacific SST anomalies seem to change the behavior of long 
waves. The spectrally filtered Hovmoller diagram is employed to examine the 
propagation and amplification of long waves. It is found that wavenumber 2 
and 3 in the anomaly run became stationary around the climatological locations 
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ol obst>rvat lotial study (Austin, 1980). The constructive Interference of 
wavtMiumber 2 and 3 in the anomaly run forms two persistent open blocking 
ridges: one is in the west coast of North America and the other in the western 

Europe . 


Since the persistent blocking ridges can be formed bv the constructive 
Interterence of stationary long waves, the next problem we have to examine is 
luHv (hose persistent blocking ridges are amplified and maintained. The spectral 
energetics Is carried out to answer this problem. Firstly, tlie total energy 
budget analysis displays that the time variation of the blocking ridge 
energetics over its Lite span Is the same as that during the typical blocking 
situation as sketched by Lejentts (1977). Spectral energetics of wavenumber 
2 and 3 shows that A^ and A.^ are supplied by A 2 ) and C(Aj,, A^) , 

respectively. Note that A^ < A 2 may be caused by C(Ar,, A^) < A 2 ). 

K 2 Is provided by ^ 2 ), while Is maintained by -CCK^, K^) . In other 

words, K,, Is furnished by barocllnic energy conversion and is supplied by 
barotropic energy conversion. 

Tung and Llndzen (1979) pointed out that the severe weather patterns of 
January 1963 and 1976-1977 winter were similar to some extent. Since the 
pattern of the pacific SST anomaly used in anomaly run Is 1976-1977 winter, 
it would be of interest to compare the spectral analysis and energetics 
between the anomaly run and these two winters. The spectral energetics 
studies of l^iin-Nlelsen (1964), Wlln-Nielsen et al. (1964), and Murakami and 
Tomatsu (1965) show that the energetics of wavenumber 2 and 3 is similar to the 
anomaly run, especially wavenumber 3 in January 1963 was vitally furnished 
by the barotropic energey conversion -C(K^, . The anomaly run began from 

January 1. We, therefore, start the spectral analysis and energetics of 
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1976-1977 winter from the same day. First of all, the spectrally filtered 
Hovmoller diagrams show that long waves (n“l-4) are stationary over most of 
this time period. The open blocking ridge located in the west coast of North 
America in this winter was mainly formed by the constructive Interference of 
wavenumber 2, 3 and 4, especially 2 and 3. The blocking ridge decayed when 
these wavenumber 3 and 4 became moving. The development and maintenance of 
wavenumber 2 and 3 in this winter Is different in some way from either in 
January 1963 or the anomaly run, and of wavenumber 2 and 3 are supported 
by C(A„, A ) and C(A , K ) respectively. On the other hand, K„ is very slgnlfi- 

/j II n n 

cantly depleted by C(Kj, k^). That is, wavenumber 3 is barotroplcally stable 
in this winter. It is opposite to the dynamical behavior of wavenumber 3 
either in January 1963 or the anomaly run. 

Hie summary provided above indicates that the necessary condition to have 
blocking ridges formed is the constructive Interference of stationary long 
waves. The other condition for the formation of blocking ridges Is how these 
long waves are amplified and maintained. The spectral energetics of blocking 
ridges shows that the amplification and maintenance of stationary long waves 
can be due to various physical processes: barocllnic and barotropic. In 

fact, it is also possible due to nonlinear interaction between cyclone-scale 
and planetary-scale waves as proposed by Gall et al. (1979) in a numerical 
experiment and observed by Hansen and Chen (1981) in diagnostic analysis of 
observational study. 

Namias (1978) suggests that the Pacific SST anomaly of 1976-1977 winter 
is an Important factor to cause the persistent blocking ridge in the west 
coast of North America. His hypothesis is confirmed by the numerical 
experiments analyzed in this study. In other words, the SST anomaly may 
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cause the statlonarlty of uitralong waves. It has been a well known problem 
in the sophisticated atmospheric model that the energy content of long waves 
are too small compared to observational study. This problem is also found 
recently in the GLAS model (Chen et al., 1981) and the operational model 
of European Center for Medium Range Weather Forecasts (Hollingsworth et al., 
1980). Although Spar (1973a, b) and Spar and Atlas (1975) show that SST 
anomaly may not improve very much the long-term forecast using the GISS 
atmospheric model, our present analysis indicates that the careful treatment 
of the SST anomaly may be important to the simulation of long waves which may 
in turn improve the long- term forecast. 

It has been observed that the stratospheric major warming always 
associated with the tropospheric blocking (e.g. Labitzke, 1967; Schoeberl, 
1979). Our present analysis and observational study of Austin (1980) show 
that the persistent open blocking ridges or Rex-type blockings are established 
due to the constructive Interference of stationary long waves. It is shown 
In the anomaly run that the SST anomalies may force long waves s tationary . 

Tlie introduction of SST anomalies to a general circulation with stratosphere 
may increase the possibility of success to simulate the stratospheric warming 
as Miyakoda et al . (1970). 
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A ppendix; 


The spectral energy equations (l)-(4) can be illustrated by the following 
energy diagram, Figure Al. The notations used have been explained In Section 
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1. Observed Pacific sea surface temperature anomalies (°F) during 
the winter of 1976-1977 (after Namlas, 1978) . 

2. TVo-week (February 1-14) mean height fields: (a) 850 mb control 

run, (b) 850 mb anomaly run, (c) 500 mb control run, and (d) 

500 mb anomaly run. 

3. Spectrally filtered Hovmdller diagrams at 500 mb and 50°N for 

three wave groups: (a) wavenumber 1-15, (b) wavenumber 1-4, 

and (c) wavenumber 5-15. 

4. Spectrally filtered Hovmbller diagrams at 500 mb and 50”N for 
four ultralong waves and combination of wavenumber 2 and 3. 

5. Spectrally filtered synoptic charts of height and temperature of 
wavenumber 2, 3 and combination of them at 500 mb on February 4. 

6. Time variations of energies: (a) and A„, (b) K„ and K^, and 

(c) energy conversions: Aj.),'C(Ag, Kj^) and C(Rp, K^) . 

7. Schematic picture of the energy changes and energy conversions 
during a blocking situation. Day t ■ t^ Is the onset of the 
blocking (after Lejenlis, 1977). 

8. Time variation of (a) zonal flow and (b) zonal temperature at 
500 mb for the anomaly run. 

9. Ti:ne variation of momejatum transport by wavenumber 2, uv(2), 
and by wavenumber 3, uv(3), at 500 mb for the anomaly run. 

10. ^le same a^ Figure 9, except for sensible heat transport, i.e. 
vf(2) and vf(3). 

11. Geographic distribution of uv, vT and mT for wavenumber 2, 3 
and combination of them at 500 mb on February 4. 

12. Time variation of (a) available potential energy and (b) kinetic 
energy for wavenumber 1, 2, 3 and 4 in the anomaly run. 

13. Spectral energetics of (a) wavenumber 2 and (b) wavenumber 3 in 
the anomaly run. 

14. Average contours of 700 mb surface (tens of feet) for January 
1963. Troughs are indicated by solid lines (after O'Connor 1963). 

15. Energetics of January 1963: (a) time variation of A^ and Ap, 

(b) time variation of K„ and K^, and (c) Lorenz energy cycle 
(after Wiin-Nielsen, 1954). Unit of energy conversions: Wm~ . 


Figure 16. Spectral energetics of January 196J: (a) Ap spectrum, (b) 

Kp spectrum, (c) CfA^, A ) spectrum, (d) C(K , K^) spectrum 
arid (e) time variation of C(A , K ). These figures are ex- 
tracted from various studies:” (a^ and (b) from Wlln-Nlelsen 
(1964), (c) and (d) from Wiin-Nielsen et al. (1964), and (e) 
from Murakami and Tomatsu (1965). 

Figure 17. Mean 700 mb contours (tens of feet) of winter 1976-77 
(December, January, and February) (after Namlas, 1978). 

Figure 18. Spectrally filtered llovmbller diagrams at 500 mb and 50“N for 
January-February 1977. 


Figure 19. Linergy cycle of January-February 1 977 . Unit of energy con- 
versions: 


Figure JO: Spectral energetics of January-February 1977: (a) Aj, spectrum, 

(b) Kp spectrum, (c) C(A„, A„) spectrum, (d) 0(A , K„") spectrum, 
(e) CTK|^, K^) specirum, (f) CA(N|m,0 and (g^ CK(n|m,)i). 


Figure Al. 


Stliematlc energy .'i.igram for wavenumber n. 
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Abstract 


The spectral energetics of two blocking case studies froa the winter of 
1978-79 are calculated. One case occurred over the North Atlantic and the 
other over the North Pacific. The tewporal evolution and geographical dlstrl 
butlon of the energetics of the planetary-scale (sonal hansonlc wavenumbers 
1-4) and the Ittermedlata-scale (wavenumbers 5-10) wsve ensembles are exam- 
ined to Identify the physical mechanisms Important In the blocking process. 

Two distinct forcing mechanisms were found. The Atlantic block was 
forced by the nonlinear Interaction of Intense barocllnlc cyclone-scale 
waves with ultralong waves. The Pacific blocking resulted from the 
barocllnlc amplification of planetary-scale waves. At least two necessary 
conditions are required for the development of blocking. First, cne of the 
forcing mechanisms must be In operation, and second, the phase of the ante- 
cedent planetary-scale waves must be such that the forcing mechanism can 
amplify them Into a blocking pattern. 

In each of the cases studied, the block developed downstream of intense 
intermediate- or large-scale cyclogenesis. The possible relationship between 
this cyclogenesis, subsequent block development, and the pattern of sea 
curface temperature anomalies In the western portions of the oceans Is 


discussed. 
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1. Introduction 

Blocking action can raault to aignlf leant local dlatortiona ol noraal 
aaaaonal climatic tronda (Bax* 1950b). Racont abnoraal racord aatting wintar« 
in North Aaarlea, particularly 1976-77 and 1978-79, have atiiiulatad 
ronawad Intaraat in tha blocking problem. Numaroua articlaa dlaeuaalng 
varlouB aapacta of thla aubject hava appeared in recent literature. A 
problem frequently encountered when atudying blocking action ia eatablishlng 
exactly what la meant by the term blocking. Any peralatent, quaaiBtaklonary 
ridge that impedea the normal eaatward progreaalon of cyclone-acale waves 
can be called a block. Our primary intereat ia in high latitude blocking 
ridgea that have lifetimea of 10 to 15 daya and that generally aatlsfy 
Rex'a (1950a) definition of blocking. 

The fact that blocking can cauae large, aeaaonal climatic anomaliea 
and that the physical mechanisms responsible for its formation are not well 
understood suggests that additional diagnostic studies of blocking are 
needed. The spectral energetics method provides a powerful diagnostic tool 
for the qualitative identification of mechanisms important in block 
development. Examination of the temporal evolution of the energetics 
can be particularly useful, especially if an effort is made to relate salient 
features of the energetics analysis to synoptically observable events. 

We have calculated the evolution of the spectral energetics of two 
cases of blocking from the winter of 1978-79. The forcing mechanisms of 
these two ridges were quite different. The first case occurred over the 
North Atlantic in December 1978. This block was forced by the nonlinear 
Interaction of barocllnlc cyclone-scale waves with planetary waves. The 
second case occurred over the North Pacific in late December 1978 and early 
January 1979. This block developed as a result of the baroclinic amplifl- 
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cation of planetary-acale wavea. Tha forcing machaniana of theaa blocking 
eplaodaa are rapraaantatlva of othar caaea of Atlantic and Pacific blocking 
(Hanaan, 1981). 

2. Background 

Early Inveetlgatcra (Berggran, at al., 1949; Elliott and Smith, 1949; 
Rex, 1950a, b) eatabllahed the aynoptlc behavior and climatology of block- 
ing. According to Rex (1950a), blocking exists when the basic westerly 
current at upper levels splits Into two branches that extend over at least 
45* of longitude, with each branch trasnportlng appreciable mass. A sharp 
transition from zonal flow upstream to meridional flow downstream of the 
current split must be observed and this pattern must persist for at least 
ten days. Rex (1950b) noted that blocking occurs roost frequently over the 
northeastern portions of the Atlantic and Pacific Oceans, normally persists 
for 12 to 16 days, and Is relatively stable In position. The climatological 
mean position of the blocks is, of course, just downstream from the normal 
position of the major midlatitude jet streams. The original single westerly 
jet Is replaced by two westerly branches with a zone of easterly flow in 
between. A pronounced warming occurs In the northern part of the blocked 
zone with cooling In the southern part. The mean meridional temperature 
gradient across the blocked zone Is greatly reduced or even reversed (Rex, 
1950b) . 

Rex (1950a) observed that the formation of a blocking anticyclone Is 
Initiated by a "finite external-impulse" provided by a sufficiently Intense 
cyclonic disturbance. Cyclone-scale waves apparently played an Important 
role In developing the Atlantic block studied by Berggren et al. (1949). 
What, If any, cause and effect relationship exists between the cyclones 
and subsequent blocking cannot be determined from synoptic studies alone. 
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lUiCcllffe and Murray (1970) found that cold aoa aurfaco tomparature 
(8ST) anomallas in a vlda araa aouth of Nawfoundland ara aaaociatad with 
blockad waathar pattarna over northarn and waatarn Europa. Warn SST 
anonallaa in tha earn araa favor nora prograaaiva flow pattarna. Naniaa 
(1964) nada a ainilar obaarvation and auggaatad that cold SST anoMliaa naar 
Nawfoundland cauaa incraaaod waatarn Atlantic cycloganaala which in turn 
laada to block davalopnant ovar northarn Europa. It ia difficult to con- 
flrn thla apaculatlon becauaa tha propoaad procaaaaa occur in data aparae 
raglona ovar the ocean. 

Recently* several new thaorlaa for blocking have bean proposer*. 

Eggar (1978* 1979) uses a barotroplc channel nodal to show that blocking 
can result from tha nonlinear interaction of forced stationary waves with 
slowly moving free waves. Charnay and Devore (1979) use a low-order, 
barotroplc channel model to suggest that for a given forcing, both high 
and low index equilibrium flows can exist. Blocking could be a 
meta-stable equilibrium of the low index type. Charnay and Strauss (1980) 
have extended this idea to a baroclinic model, where they found that oro- 
graphy played an indirect, catalytic role in developmsnt of a low-index, 
blocked state. The transition to this blocked state may occur, in their 
opinion, through a combination oi orographic and baroclinic instability. 

Tung and LlAdaan 0.979) show analytically that blocking may be caused by 
linear, resonant amplification of planetary-scale waves forced by topog- 
raphy and surface heating. If the atmospheric flow is such that it renders 
a wave stationary with respect to the earth's surface, the wave can inter- 
act resonantly with the stationary forcing. 

None of the recent theories specifically Incorporate the forcing of 
planetary-scale waves by cyclone-scale waves. However, Gall et al. (1979) 
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hav« BugiMtad that the nonuiilfor* •pacing and intanalCy of cyclonaa around 
a lacltuda sona will lead to forcing of the ultralong wavea by axially 
aaynBMitric cyclone-acale heat tranaporta. The reaultlng pattern of high 
latitude warming and lower latitude cooling ahown by Gall et al. reaeiri>lea 
that aaaoclated with Rex blocking. The aaywietrlc cyclone-acale heat 
tranaporta In turn force large-acale vertical motion and correaponding 
large-acale barocllnlc energy converalona that provide the kinetic energy 
source for ultralong wave development. 

Previous energetics analyses of blocking, using primarily time-averaged 
results, indicate that blocking la accompanied by strong barocllnlc energy 
conversions. Julian and Labltske (1965) noted Increased vigor In the 
vertical motion field during the onset of blocking and concluded that the 
energy variations associated with It were largely barocllnlc. Paulin (1970) 
and Murakami and Tomatsu (1965) found that large conversions of sonal to 
eddy available potential energy and, then to eddy kinetic energy accom- 
panied block development. The maximum monthly mean barocllnlc conversions 
occurred at zonal harmonic wavenumbers 2 and 3. Winston and Krueger (1961) 
discovered similar barocllnlc energy conversions during the development 
phase of an Atlantic block. 

All the aforementioned studies use data from the entire middle and high 
latitudes of the Northern Hemisphere. Zonally averaged results like these 
cannot distinguish whether diagnosed processes are associated with blocking 
or with some other aspect of the circulation. In addition, the calculation 
of monthly means can easily mask Important energetics processes because the 
lifetime of blocking is usually only 10 to 15 days. 
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Nor« rsemtly* Hurtiiaii and (ftian (1980) uaad a atatiatical approach to 
axaMlna tha haat and vorticlty btidgata of conpoalttd Atlantic and Pacific 
blocka. Thay auggaatad that Pacific blocka aay ba cauaad by largaly baro- 
tropic aMchanlaMt whila Atlantic blocka my ba largaly baroclinic. Auatln 
(1980) haa ahown that blocking can ba rapraaantad by tha conatructlva 
inter faranca of planatary-acala xonal harmonic wavaa with normal phaaaa 
but much graatar than normal amplitudaa. 

Howavar* ao far no machaniam haa baan iaolatad to explain tha initial 
development of blocking. Although it appaara from moot atudlao that baro- 
clinic procaaaea are important* tha temporal evolution fo tha apectral 
anergatica of blocking haa not yet baan prarantad. 

3. Data and Procedurea 

Tha aquationa governing atmoapharic anergatica in wavenumber domain 
were firat formulated by Saltzman (1957). The complete apectral energetics 
formulation can be found in numaroua other aourcea (e.g., Hanaen, 1981; 
Saltzman, 1970) and will not be repeated here. Note that the nonlinear 
interaction represents a redistribution of energy between the various waves 
and is neither a source nor a sink of kinetle. energy. Whan summed over all 
wavenumbers* tha nonlinear interaction should yield no net change in total 
eddy kinetic energy (Steinberg et al.* 1971). We have reformulated the 
spectral energetics equations in such a way that the wavenumber domain equa- 
tions can be reconciled with the physical domain equations* and the nonlinear 
Interaction terms sum to zero. Our reformulation of the nonlinear inter- 
action terms is given in the Appendix. 
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Th« data uaad In this study ara tlia output of tits National Notaoro- 
logical Cantar'a (NMC) operational analysis. Twice daily Northern He«i- 
sphara data for the horlsontal wind (u, v), the gaopotantial height (*), 
and the tanperatura (T) on a 2.5* by 2.5* grid at tha 10 mandatory levels 
In the troposphere (1000, 850, 700, 500, 400, 3C3, 250, 150 and 100 mb) 
were obtained from the National Center for Atmospheric Research. In 
order to measure the release of available potential energy (APE) i trough 
baroclinlc energy conversion, it is necessary to compute the vertical 
motion in the atmosphere. The vertical velocity in pressure coordinates 
(w) was calculated from a spectral form of the quasigeostrophlc w-equation 
assuming w - 0 at 1000 mb and 100 mb (Chen et al., 1981). Fourier 
coefficients of u, v, T and w at each latitude circle and every level 
were calculated using a standard fast Fourier transform technique. The 
zonal harmonic expansion of the dynamical variables was truncated after 
wavenumber 18. 

The spectral energetics results were Integrated over the latitude 
zone from 30*N to 80*N. This zone was selected because It completely 
Includes the features of the middle and high latitude circulation associ- 
ated with blocking, while excluding thv Influences of the tropical Hadley 
cell circulation. The meridionally averaged results are also Integrated 
in the vertical from 1000 to 100 mb. 

In order to more easily relate the spectral energetics results to 
observable features of the circulation, ensembles of planetary-scale (wave- 
numbers 1-4) and Intermediate-scale waves (wavenumbers 5-10) were formed. 

A similar approach has been used in previous energetics studies (Wiln-Nlelsen 
ot al., 1963, 1964; Chen et al. 1981). By examining the results for these 
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•naemblcii, %m can laolate acnoapharlc cncrgatics proceaaaa occurring on 
planetary or intarnadlarc acalea an well aa nonltnaar Interactiona between 
theae acalea. Sauiller acale proceaaea (wavenuabera jn 11) were found to play 
no aignlf leant role In blocking developaent. 

It la not poauible to detemlne froa the tlae evolution of the apectral 
energetlca alone whether the dlagnoaed energy or energy converalona have 
anything to do with the developaent of blocking. He ahould not expect that 
the developaent or exlatence of a blocking high will neceaaarlly have a 
profound effect on the sonally averaged apectral energetlca. Contrlbutlona 
from other regions can easily mask the energy processes related tc blocking. 
We need to determine the spatial dlatrlbutlon of the energetlca results In 
order to draw more definite conclusions. 

Polar stereographic projections of planetary- and Intermediate-scale 
ensemble kinetic energy, transports, height and temperature fields and 
barocllnlc conversions were made for Individual observation times when 
energy values or energy conversions were maximum or rapidly changing. 

These plots provide considerable Insight Into the distribution of energetics 
processes around the hemisphere and make It much easier to identify impor- 
tant processes related to blocking. To supplement the spectral energetics 
calculation, limited area energetics calculations were also performed. 
Specific Information we wish to obtain from the limited area calculation 
Is the location of wave-wave energy transfers. 

4. A Case of High Latitude Atlantic Blocking 
A. Synoptic Behavior 

During the second half of December 1978, a major blocking ridge rapidly 
developed over the North Atlantic Ocean. Preceding the development, strong 
westerly flow between 40"N and 50*N dominated the Northern Hemisphere 
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circulation. A ohort wave 500 mb trough appeared over the weatern Great 
Lakea on 17 December In a region of strong aonal flow. This feature 
developed rapidly Into a significant northwest-southeast tilted trough 
with strong warm advectlon apparent on its downstream side. By 0000 GMT 19 
December, the trough had deepened considerably into a 494 dam closed low 
centered Just north of Maine, with a trough extending southeastwa :d over the 
Atlantic. A ridge was building south of Greenland ahead of this trough. 

The trough had acquired a very pronounced northwest-southeast tilt by 
0000 GMT 21 December (Fig. 1) as the ridge continued to develop. 

The block reached maturity at 0000 GMT 22 December as a 568 dam closed 
high centered at 68”N over the Greenland shore of the Denmark Strait. The 
original trough was still northwest-southeast tilted but It had greatly 
weakened. A stable, blocking high remained over Greenland centered near 
70"N with 500 mb heights ranging from 560 to 569 dam until the 28th. Strong 
500 mb easterlies along 60*8 persisted throughout this period. At 1200 GMT of 
the 25th, easterlies extended from the Norwegian coast to Baffin Island and 
the 500 mb pattern had a classic Rex blocking structure (Fig. 2). 

Beginning on roughly 28 December, high latitude retrogression of the 
blocking ridge was evident. At 1200 GMT 29 December (Fig. 3) the 500 mb 
pattern still resembled Rex blocking, but extension of the 500 mb ridge over 
Baffin island and the Queen Elizabeth islands was apparently the result of 
retrogression. Further height rises In this vicinity, particularly over 
the Beaufort Sea occurred between the 29th and the end of the month. By 
0000 GMT 1 January 1979, the Atlantic block had completely decayed. 

The development phase of the block was from roughly 19 to 22 December, 
rho mature phase from 22 to 29 December, and decay occurred rapidly from 
29 December to 1 January. The synoptic behavior of the case is very 
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similar Co Chat described by Berggran ec al. (1949). A speccrally filtered 
HovmoUar diagram of che plane Cary-scale 500 mb height field at 70*N 

clearly shows Che' strong ultralong wave ridge amplifying during 19-22 
December and persisting until the end of the month (Fig. 4) with retro- 
gression apparent after roughly 27 December. The mean longitude of the 
high center is about 30*W. As we shall see, the largest energy conversions 
in Che entire hemisphere at this time were associated with the development 
of this block. 

B. Energetics 

Before examining the various aspects of the ensemble results for this 
case, let us first consider the time evolution of the overall energy cycle. 
The vertical Integration of the zonal kinetic energy (K„), eddy kinetic 
energy (Kg), zonal APE (A^,) and eddy APE (Ag) averaged between 30*N and 80*N 
are shown In Fig. 5a. The conversions between these quantities are given In 
Fig. 5b. 

Kg had maintained a relatively large, steady value until 14 December 
corresponding to the aforementioned strong zonal flow In mid-December. 
dropped sharply from the 14th to the lAth, eventually reaching a minimum on 
the 20th. Minimum values of C(Ag, Ag) and C(Ag, Kg) on 14 and 15 December 
correspond to a period of rapidly Increasing A . However, beginning on 
15 December, both C(Ag, Ap) and C(Ag, K^) Increased dramatically, reaching 
striking maxima on the 17th and 18th, respectively. They remained relatively 
large until the 22nd after which they declined sharply. As a result of 
these barocllnlc processes, Ag declined rapidly, and Kp and Ag rose from 
minima on 15 and 16 December to marked maxima on the 19th and 21st, 
respectively. Notice the negative correlation between Ag and Kg (Fig. 5a) 
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that Is consistent with these conversions. Lejen&'s (1977> has noted a 
similar energy cycle during blocking In the work of previous Investigators. 

The decline In Kj,, rapid growth and decay of Aj, an4 mlnlsia in Kp and 
Ag precede the onset of blocking. Very large barocllnlc processes 
(C(Ap, Ag) and Kg) ] just prior to and during the development stage 

of the block account for the observed Increases In Ag and Kg In the third 
w* t'k of December , with a much smaller contribution to K„ from C(K„, K„) . 

E Cl 

The largest eddy kinetic energy values correspond to the developing phase 
of the block. Kg declines and remains relatively low during the mature 
stage. Ag reached its maximum as the block reached maturity. 

The ensemble results give a clearer picture of the energetics of this 
case. The vertically integrated energetics of the planetary-scale ensemble 
(wavenumbers 1-4) and the Intermediate ensemble (wavenumbers 5-10) are 
shown In Figs. 6a and 6b. 

The large values of C(Ag, and C(Aj_^ *^1-4^ 

December and the CfKg, Kj^_^) maximum on the 18th (Fig. 6a) are primarily 
associated with a deepening ultralong-wave trough over East Asia. How- 
ever, the intermediate-scale kinetic energy and energy conversions from 
I'i to 22 December (Fig. 6b) were concentrated over eastern North America 
and the western Atlantic. The intermediate-scale kinetic energy (K^__j^q) 
increased dramatically from a minimum on 14 December to a striking 
maximum on the I9th. accounted for roughly half of the total eddy 

kinetic energy which was also maximum on that date. The Intermediate- 
scale .\PE (A 5 _|q) underwent a similar but less dramatic increase. These 


energy Increases are the result of large barocllnlc energy conversions. 


The larger magnitude of CCAp. 


A5-10) than K5 _^q) 


accounts for the 
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increase In A« from 16 to 18 Deceiriier. Note that the tn termed la te- 
5-10 

acale converalona repreaent moat of the total Aj, to Ag and Ag to Kg con- 
versions before and during the block's development. The decline in Ag 
from 16 to 19 December has a marked negative correlation with the increases 
in Aj_j^jj and during the same period. 

The major features of the intermediate-scale 500 mb height field at 
1200 GMT 18 December (Fig. 7a) were a deep, northwest -southeast tilted trough 
extending from eastern North America into the western Atlantic with strong 
cyclone-scale ridges both upstream and downstream of the trough. The 
position of the trough over eastern North America in the synoptic chart 
for this date corresponds exactly to the intermediate-scale trough position 
in Fig. 7a. The cyclone-scale barocllnlc energy conversion was maximum 
at this time and C(Aj_j^q, ^ 5 .^ 0 ^ concentrated around the eastern North 
American cyclone system (Fig. 7b). The largest conversion is over 
Laborador in southeasterly flow on the downstream side of the trough. 

C(Ag, Aj_jq) was also very large at 1200 GlfT 18 December. The cyclone- 
scale heat transports and therefore C(Ag, were concentrated arotmd 

the same trough-ridge system. Northward transport of warm air east of the 
trough axis and southward transport of cold all «^est of the trough axis 
can be inferred. By 1200 GMT 19 December, the heat transport associated with 
this trough had Increased, especially the northward warm air transport on 
the trough's east side, and the convergence of the cyclone-scale heat trans- 
port had reached the south tip of Greenland. The zonal asymmetry of the 
cyclone-scale heat transport was particularly striking at 850 mb (Fig. 8). 

The northwest-southeast tilting of the trough-ridge system caused con- 
siderable southward transport of westerly momentum (northward transport of 
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eascerly momfititum) . The negative cyclone-scale nomentum transport at 
1200 GMT of the 19th was especially promlnant over Newfoundland and south of 
Greenland (Fig. 9a). The ultralong wave 500 mb height field at this time 
(Fig. 9b) shows a diffuse ridge over the North Atlantic with the suggestion 
of easterly (geostrophlc) flow from south of Greenland to over Hudson Bay 
(a similar pattern existed at 300 mb). Thus from 18 to 22 December there 
was convergence of cyclone-scale easterly momentum Into a region of 
planetary-scale easterly flow (compare Figs. 9a and b). This coincides wltli 
the period of large nonlinear interaction between the Intermediate and 
planetary-scale waves. A clear negative correlation between the planetary- 
scale and Intermediate-scale kinetic energy existed from 19 to 23 December, 
consistent with the diagnosed nonlinear Interaction (compare Figs. 6a and b). 
Tlie period of declining and large nonlinear transfer of energy to the 

rapidly rising from 19 through 22 D'icember corresponds exactly to the 

time of the rapid development of the North Atlantic blocking. 

Comparison of the Intermediate-scale and planetary-scale kinetic energy 
distributions over eastern North America and the North Atlantic on the 19th 
and 23rd supports the contention that Intermediate-scale kinetic energy was 
transferred nonllnearly to the ultralong waves In this region. On the 19th, 
*^5-l0 ^ very remarkable, zonally assanmetrlc distribution In this area 

(Fig. lOa) but was nearly nonexistent. By 1200 GMT 23 December, this 

pattern had reversed. The maxima over the Atlantic had completely 

disappeared. There are no contours over eastern North America or the Atlantic 
In the Plot at the time (not shown). Howevei , very large values, 

associated with the blocking high over Greenland, had replaced the inter- 
mediate-scale kinetic energy (Fig. 10b). Easterly flow south of Greenland 
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was particularly strong. Thus, ths mechanism for the large diagnosed 
energy transfer from intermediate to planetary scale was apparently the 
northward transport of easterly momentum by cyclone waves into a region of 
planetary-scale easterly flow over the North Atlantic, resulting In the 
acceleration of the easterly flow. 

At 1200 GMT 19 December, when the blocking ridge was rapidly developing, 
«*!-«• *“• «>. m,.t part. the North Pacific 

sector and a region just south of Greenland (Fig. 11). The large-scale ascent 
of warm air south of Greenland could have been forcod by the strong cyclone- 
scale heat flux convergence In this area in a manner similar to that pro- 
posed by Gall et al. (1979). In any event, the source of provided 
by the nonlinear Interaction was as large or larger than the harocllnlc 
source when the block was developing from 19 through 22 Deceiriber (Fig. ha), 
but the baroclinic conversion was distributed in various locations during 
this time. Therefore, we can conclude that the nonlinear interaction was 
the primary source for the kinetic energy of the block. 

The equivalent barotroplc structure of the ultralong waves over the 
Atlantic during the block's development is shown by comparing the planetary- 
scale height and temperature fields. At 1200 GMT 19 December , the weak 
ridges in the 500 mb height and temperature fields over the Atlantic were 
in phase. Even after the block had reached maturity at 1200 GMT 23 December, 
the amplified and Tj^_^ ridge axes over Greenland were still almost 
exactly in phase at 500 mb (compare Figs. 12a and b), as well as at 850 mb 
and 300 mb. Notice the very large Increase in the amplitude of the ridge 
over Greenland from the 19th to the 23rd (compare Figs. 9b and 12a). The 
thermal ridge experienced a correspondingly large amplification, consistent 
with the increase in 
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A maxi»um In the planetury-scale APE on 22 and 23 Oecaabar corraaponds 
to the time when the thermal ridge aaaociated with the Atlantic block had 
Ita greatest amplitude. A major source of this APE was provided by the 
excess of C(Ajji over C(Aj^__^» between 1200 GMT 19 December and 

0000 GMT 23 December (Fig. 6a). The CCAj^t maximum on the 22nd was 

associated with significant net positive planetary-scale heat transports In 
the vicinity of the block (Fig. 13), although part of the marked northward 
warm air transport centered over Baffin Bay at 500 mb was counterbalanced 
by southward transport Just east of Greenland. This viewpoint Is supported 
by the build-up, noted earlier, of the planetary-scale thermal ridge over 
Greenland between 19 and 23 December. 

An additional source of Aj,_^ during the development phase was the non- 
linear transfer of intermediate-scale APE to the ultralong waves. Normally, 
the planetary waves export APE to smaller scales but a small positive con- 
tribution of to was maintained from 9 December to 21 

December (not showiO • The cyclone-scale heat transports over the western 
Atlantic with their attendant Aj, to A^_jq conversion were apparently 
followed by a transfer of some A^_jq to the planetary waves. (Note the 
decline in A^_j^q during this time in Fig. 6b). The positive to Aj^ ^ 

transfer in this region was evidently large enough to more than counter- 
balance the normal negative transfers that were most likely occurring in 
the rest of the hemisphere, thus yielding the net positive transfer in the 
zonal ly averaged energetics. 

The ni'nlineat interactions that forced the development of the Atlantic 
block dropped to zero on 23 December by which time the planetary-scale 
kin«-Llc energy was very large (Fig. 10b) and an amplified ultralong wave 
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rtdg^ was In place over Greenland (Fig. 12a). Notice the very obvious 
wavenunber 1 pattern north of 60*N in Fig. 12a. 

All of the large-'scale energy conversions fell to nearly zero on 23 
and 24 Decenber and declined rapidly (Fig. 6a). No significant net 
energy sources are apparent in the vicinity of the block throughout Its 
mature phase as the kinetic energy around the block declined markedly. 
However, the amplitudes of the planetary-scale height and thermal ridges 
over Greenland remained high from 22 to 28 December. The and 
ridges remained approximately in phase during this period. 

By 1200 GMT 28 December, the block had begun to retrograde across 
northern Canada. The Greenland thermal ridge was still fairly strong on 
the 28th but by 30 December the effects of retrogression had greatly 
diminished the Greenland block. Southward transport uf warm air by the 
large-scale flow at 850 mb and 500 mb Just south of Greenland on 30 
December may also have aided the decline of the block. This southward lu.it 
tronsport was the result of the retrograde height wave moving out of phase 
with the thermal wave. 

In the last few days of December, the large planetary-scale energy 
conversions (Fig. 6a) were associated with developments in the Pacific sec- 
tor that will be discussed In the next section. Intermediate-scale pro- 
cesses were unimportant during the block's mature and decay phases. 

Corroboration of the nonlinear forcing of the block's development can 
be gained from the limited area energetics. We calculated the energetics 
within the sector bounded by 85"W, 20"W, 35"N and 75"N. This region contains 
the blocking high as well as the cyclone system that preceded it. The 
maxima In Kg and Ag on the I9th, and large values of C<Ag, Ag) and 
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C(Ag, Kg) from 18 to 22 December within the limited .tren correepond very 
closely to the Intermedlate^scele ensemble results (compare Fig. 14 with 
Fig, 6b). The large east-west flux of Kg around the I9th was largely the 
result of cyclone-scale Kg entering the western boundary of the limited 
area. The very large transfer of eddy kinetic energy to the mean flow that 
folUtwed the Kg maximum corresponds to the nonltne.ir transfer of Interraediate- 
sialo kinetic energy to the planetary waves In the spectral calculation (again, 
t’omrare Fig. 14 to Fig. 6b). The area used In the limited area energetics 
calculation equals roughly 1/6 of the spectral domain, and the magnitude 
of C(K„, K„) within the limited area was about six times the magnitude of 
the nonlinear Kg transfer in the spectral calculation. Thus, the limited 
area energetics confirms that t l»o location of the large nonlinear Interac- 
tion from 18 to 23 December 1978 was over eastern North America and the 
western Atlantic. This result strengthens our argument that the resulting 
blocking ridge was forced primarily by this nonlinear interaction. 

In summary, the development of this case of Atlantic blocking was 
preceded by a high index circulation with a large baroclinlc instability as 
indicated by a large value of A^. The nonlinear transfer of kinetic energy 
from barocliv'.tcally active cyclone-scale waves to the planetary-scale waves 
provided the major source of kinetic energy for the block. The marked 
conat asymmetry of the cyclone-scale heat transports coupled with large- 
scale zonal to eddy APE conversion caused the amplification of the thermal 
ridge associated with the block. During its mature stage, no active energy 
conversions were apparent in the vicinity of the block and Kj^_^ declined. 

The decay of the Mock was directly linked to the westward retrogression 
of the height wave. Slight baroclinlc damping may have aided the block's 


licnifse. 
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The energy cycle during the developaenc of rhle nonllnearly forced caee 
can be BUMMrlaed achematlcally aa followat 



The dotted arrowa represent relatively smaller magnitude energy conversions. 

5. A Case of Pacific Blocking 
A. Synoptic Behavior 

Development of a major blocking ridge over Alaska and the eastern 
Pacific In late December 1978 and early January 1979 closely followed the 
development of the December 1978 Atlantic block. The Pacific block <!ould 
be viewed as a westward redevelopment of the earlier feature, but we will 
treat It as a new development, keeping in mind the possible ImportantH* of 
the retrogression of the blocking high across northern Canada. 

On 25 December 1978 westerly flow at 500 mb existed between 30*N and 
60*N from central Asia to south of Kamchatka where the westerly current 
narrowed. The strongest westerlies in the eastern Pacific were between 
50*N and 60*N across the Bering Sea and the Alaskan Gulf. A mean ridge 
was present over the eastern Pacific off the U.S. West coast. By 29 
December a ridge had amplified along 140*W extending from 45**N northward 
over Che Alaskan peninsula, and an ultralong wave trough was developing over 
East Asia. A major closed low had fonned in southwestern Canada in response 
to cold air plunging southward on Che lee side of Che Pacific ridge. At 
1200 GMT 30 December (Fig. 15), there was further indication of regrogres- 
slon of high 500 mb heights across northern Canada and the large-scale 
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trough ovor eascern Siberia had deepened noticeably. The pattern over the 
Pacific quallflee as blocking on this date. 

A closed 557 dam high was positioned over Alaska at 0000 1 January 

1979. In the previous three days, the East Asian trough had progressed 
<'ii:itward out of Siberia and was centered Just east of Kamchatka on the first. 

A closed 557 dam high was positioned over Alaska at OOE 1 January 1979. 
In the previous three days, the East Asian trough had progressed eastward 
out of Siberia and was centered Just east of Kamchatka on the first. 

Between t and 4 January further amplification of the Alaskan high took 
place. The pattern during this period can be called an omega block because 
of the de»'p ultra long wave lows over the Bering Sea and Hudson Bay flanking 
the high centered at 65®N, 140*W (Fig. 16). 

The 500 mb helgitt over Alaska remained at or above 560 dam from 4 to 8 
January, llte ridge developed a northwest-southeast tilt during this time and 
movement of the Bering Sea low to south of the Aleutian Islands gave the 
p.ittern the appearance of strong, Rex blocking. 

Another ultralong wave trough began to progress slowly eastward trom 
northeast Siberia on or before 8 January, reaching a maximum depth of 472 
dam just off the Slbei.'f.n coast on 11 January. Concurrently, the block 
began to decay and by 12 January it has weakened considerably (Fig. 17). 
Westerlies bad reestablished themselves across the Pacific and the high 
center had moved north of Alaska and lost amplitude. Remnants of the ridge 
completely disappeared over the next 3 to 4 days. 

During the first week of January when the Pacific block was in its 
mature stage, strong cyclogenesis occurred over eastern North America 
tul lowed by the developiuenl of a major midlatitude ridge over the Atlantic. 


Hie longitude vereiie time eectlon of the planetary-scale height field 
(wavenuMbera 1-4) at 60*N (Fig. 18) shows the stationary poaitloii of the 
blocking ridge along 140*U througliout its life cycle. Note the troughs 
both upetreasi and downstream of the blocking ridge during Its mature stage 
<30 Dec. to 8 Jan.). 

B. Energetics 

The time history of the vertically integrated energy cycle for the 
sonal mean and eddy components of the flow during the December 1978-January 
1979 Pacific blocking are given in Figs. 19a and b. Steadily increasing 
barociinic processes consistent with rising Kg and declining Ag, coupled 
with a relatively large C(Kg, Kg) nark the energy cycle during the block's 
developing phase In the last few days of December. C(Ag, Ag) and C(Ag, Kg) 
Increased steadily from 25 December to 1 January, while C(Kg, Kg) was 
strongly positive from the 28th to the 1st. 

Moat of the dramatic energy conversions that occurred In the zonallv 
averaged energy cycle from 2 through 8 January were related to developments 
downstream from the Pacific ridge. 

After 8 January as the Pacific block decayed. Kg rose slowly while Ag 
declined, "nie negative value of C(Kg, Kg) during this period may explain 
the Kg increase. Only a slight increase in Ag and Kg accompany the 
C(Ag, Ag) and C(Ag, Kg) maxima on the 9th. Notice that Ag is negatively 
correlated with both Ag and Kg throughout the block's history. This rela- 
tionship emphasizes the importance of barociinic processes in the zonally 
averaged energetics during this period. 

The planetary waves account for the large majority of the energy total 
conversions shown in late December. These energy conversions were primarily 
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assoclaCwd vlch davclopMnti over th« North Pacific and Ea^t Aala. Tha 
planetary-'Mcalc kinetic energy waa ateady from 26 through 23 Oecetiiber (Fig. 
2U4)» but it roae aharply from the 29th to the 31at. Large-acale baroclinlc 
proceaaea IC(Ajj, A^_^) and C(Aj^_^» were very large in the laat aix 

daya of December. On the 26th and 27th* negative CCK^. and wave-wave 

Interactlona countet balance the baroclinlc kinetic energy generation. How- 
ever, from the 28th through the 31at, CfK^. became atrongly poaitlve 

unci the nonlinear interaction remained near aero allowing the very large 
C(Aj_^, and the large to cauae the large-acale kinetic 

energy to Increaae aharply. 

Subatantial southward transport of cold air under the ultralong wave 
East Asian trough and northward transport of warm air over Alaska and the 
Bering Strait had developed by 28 December . The horizontal plots for 1200 
CMT 30 December are representative of the 28 to 31 December period when 
large-scale baroclinlc and barotroplc processes were quite large. The 
East Asian planetary-scale low was centered just south of Kamchatka at 1200 
GMl' 30 December and the retrograde ridge over Alaska had begun to amplify 
(Fig. 21a). C(Ajj, was near a marked maximum on this date (Fig. 20a) 

and the most significant ultralong wave heat transports through the depth 
of the troposphere were over the Bering Sea and along the Aslan coast. The 
northward transport of warm air Is particularly strong at 850 rob and 500 mb 
over the Bering Sea (Fig. 21b). Comparison of the heat transport with the 
500 mb planetary-scale ensemble temperature pattern over the Pacific (Fig. 
21c) shows that the convergence of the planetary-scale eddy heat flux is 
acting to furcher amplify the thermal pattern by transporting warm (cold) 
air into the region of the thermal high (low). This is consistent with the 
large observed ClA^, 
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Th« larg* •»€«!• barocllnlc anargy convaraion at 300 «b (Fig. 2 Id) sliowa 
a vary atrong anargy convaraion in aacending vara air over the Baring Sea. 

A broad area of poaltlve energy conversion la aaaoclated with sinking cold 
air beneath the Esst Aslan trough. Comparing the 500 mb planetary-scale 
height and temperature fields at 1200 GHT 30 December (Flga. 21a and 21c) * we 
sea that the phase of the temperature field Is shifted to the %rest of the 
height field. (The (^hase shift Is more obvious in the lower troposphere 
(not shown).) This phase shift Is consistent with the observed heat trans- 
ports and C(Aj^_^f ^1-4^ noted above, and clearly demonstrates the barocllnlc 
structure of the ultralong waves over the Pacific during the Initial develop- 
ment of the block. 

In addition, the planetary-scale ensemble's 300 mb momentum transport 
(Fig. 22) shows strong northward transport between 25"N and 40“N souti* and 
east of Japan on the south side of the East Asian trough. This suggests 
that the majoi ry oS the diagnosed transfer of to the ultralong waves 

occurred In cor ; auction with the Intensification of the East Asian lew. 

The effect of these energy conversions on the planetary-scale kinetic 
energy budget Is shown In the distribution over the Pacific at 1200 GMT 
31 December (Pig. 23) when was at a maximum. The large values of 
in the vicinity of the blocking high over the Aleutian Islands, the Bering 
Straits, and north of the Bering Straits are priiarily from the barocllnlc 
sources. The very large values In east of Japan associated with inten- 
sification of the Japan Jet are due to both barotroplc and barocllnlc pro- 
cesses. The corresponding distribution on 28 December (not shown) had 
only very weak maxima over extreme northeastern Siberia and southwest of 
Japan. Intermediate-scale processes (Fig. 20b) were of no discernable 
Importance In the Pacific sector. 
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Further amplification of the block occurred between 31 December 1978 
and 4 January 1979. ^1-4^ remained relatively strong during this 

period, reaching another marked maximum on 3 January, but 
and CCK^i dropped sharply on 31 December and 1 January, respectively 

(Fig. 20a). Very strong planetary-scale heat transports existed over 
Alaska and the Bering Straits during this period on the upstream side of 
the blocking ridge (Fig. 24). Warm air was transported northward Into the 
region of the thermal high associated with the block (Fig. 25a). The 
greater amplitude of the Alaska thermal ridge at 500 mb on 4 January (Fig. 
25b) compared with 30 December (Fig. 21c) shows the cumulative effect of 
the large-scale heat transports and the resultant positive ^(Aj,, A^_^) In 
this region. Comparison of the planetary-scale height and temperait^^e 
fields at OOZ 4 January (Figs. 25a and b) shows that there Is still a slight 
phase shift of 10-15" longitude between the two. The phase shift Is more 
noticeable at 850 mb (not shown). The continued existence of the phase 
lug is consistent with the planetary-scale heat transports, and Is further 
evidence of the baroclinlc nature of this blocking ridge. 

Nothing particularly unusual had occurred in the intermediate-scale 

energy budget until the significant S-10 maximum on 3 January (Fig. 20b). 

This peak was caused by baroclinlc processes and was followed by a very 

large transfer of intermediate-scale kinetic energy to the planetary waves. 

The large buildup of cyclone-scale kinetic energy in early 

January took place over eastern North America and the Atlantic. A large 

Increase in K, , over the Atlantic sector (not shown) Is evidence of this 
1-4 

interaction. Thus, the large, diagnosed nonlinear Interaction during the 
hlook's mature phase had nothing to do with maintaining the Pacific block. 


Th« sharp drop in and Incrsases in and froa 2 to 8 January (Fig. 
21a) vara asaociatad with these Atlantic davelopaenta. 

Meanwhile, all large-scale energy conversions and transfers fell to 
near sero on 6 January as dropped sharply to a nininun on the 9th. 
Offsetting positive and negative heat transports (and C(A^_^» on the 

upstream and downstream sides of the high were typical during the block's 
mature stage. However, by 7 January southward warm air transport over 
northwestern Canada exceeded the positive contribution over the Bering 
Straits. The axis of ridge in the height field over Alaska had moved 
slightly to the west of the thermal ridge axis resulting In the observed 
net southward heat transport. The amplitude of the large-scale height and 
temperature waves on the 7th was quite large, but by 9 January, the maxi- 
mum height (b^_^) near Alaska had moved north of the peninsula and had 
fallen by about 100 meters. 

The large values of and C(Aj^_^ on 8 through 10 

January were associated with the developing, progressive ultralong wave 
trough over eastern Siberia and the western Pacific. This trough develop- 
ment caused the transition of the large-scale Pacific pattern from strong 
blocking on the 8th to a pattern dominated by an extensive, northwest- 
southeast tilted East Aslan trough on the 13th. 

In summary, the initial development of the Pacific block in late 
December 1978 was the result of very significant, planetary-scale baro- 
clinlc processes CC(A 2 » **'d C(Aj^_^ the 

development was triggered by the retrograde surge of high heights across 
northern Canada from the North Atlantic in late December. Additional 
amplification of the block from 31 December to 4 or 5 January 1979 was due 
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to large-scale heat transports. During the block's mature phase (4 to 8 
January), the zonally averaged energy conversions are quite small and no 
active energetics processes appeared to be occurring in the vicinity of 
the block. Decay of the block may have been partially due to barocllnic 
and barotroplc damping. The baroclinlc development of a prograde East 
Asian ultralong wave trough coincided with the block's decay phase. This 
active trough development and its eastward progression in particular appears 
to have destroyed the energetically inactive block. 

The energy cycle during the development of this case of Pacific Mock- 
ing can be summarized schematically as follows: 



Here again, the dotted arrow represents a smaller magnitude conversion, 
h. IHscusslon 

There have been a number of suggeslions In recent literature that 
intense cyclogenesis ciin act to force ultralong wave development. Earlier 
observational studies have also indicated an apparent relationship between 
intense cyclogenesis and blocking (Berggien et al., 1949; Rex, 1950a). 

Gall et al. (1979) have shown in model simulations that zonally asymmetric 
cyclone-scale heat transports can amplify planetary-scale waves through a 
type of nonlinear interaction. Namias (1980), Harnack (1980), and Sanders 
and Gyakum (1980) have all noted a correlation between intense cyclogenesis 
in the western and central North Atlantic and the recurrent strong blocking 
over the Greenland-lceland area that occurred during the winter of 1978- 
79 . We have found that the blocking high over Greenland in December 1978 
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vatt In fact* forced by Intenae, barucilnic cyclone-scale waves. A block 
over Greenland In Mid to late January 1979 was forced In a very slnllar 
manner (Hansen* 1981). 

The very Intense cyclogenesis over the western and central North 
Atlantic during the 1978-79 winter was enhanced by the presence of very 
strong sea surface temperature (SST) gradients off the east coast of 
North America (Sanders and Gyakum* 1980) that resulted from significant 
SST anomalies east of Newfoundland and Nova Scotia (Harnack, 1980). 
Increased barocllnlclty in the ocean automatically induces atmospheric 
baroclinlclty at the surface and through the depth of the lower atmosphere 
that Is in convective exchange with the surface (Bjerknes, 1962). A 
larger than normal SST gradient can lead to very intense cyclogenesis In 
the vicinity of the largest gradient (Sandeto and Gyakum, 1980). To 
maximize the SST gradient and therefore the low level barocllnlc Instabil- 
ity, the SST anomalies must be in phase with the normal SST distribution. 

The probability that blocking can form over the North Atlantic is 
Increased If cold SST anomalies exist east and south of Newfoundland 
(Ratcllffe and Murray, 1970). The position and strength of the cold SST 
anomalies In December and January given by Ratcllffe and Murray, based on 
several years with significant negative anomalies, are very similar to the 
cold SST anomaly in the same area for the 1978-79 winter given by Harnack 
(1980). In fact, the 1978-79 anomaly appears to be colder than Ratcllffe 
and Murray's mean. These SST anomalies create a larger than normal SST 
gradient off the east coast of North America which in turn leads to an 
enhanced barocllnlc rone in the atmosphere in this area. Explosive 
cyclogenesis with very strong northward warm air flow on the east side 
of the trough is spatmed by this baroclinlclty. A blocking anticyclone 
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develops north of the cyclone wave due to the Intense » conally aeymietrlc 
cyclone-scale heat transports and nonlinear transfer of kinetic energy from 
intermediate to large-scale waves. Thus, the explosive cyclogenesis that 
results from the enhanced baroclinic roue in turn leads to stronger cyclone- 
scale forcing of the planetary-scale waves and subsequent block development. 
Namias (196A) has suggested that a similar sequence of events can lead to 
blocking. 

However, strong cyclogenesis followed be substantial nonlinear 
transfer of kinetic energy from Intermediate to large scales does not Insure 
the devek'pment of blocking. For example, the development of the short dura- 
tion (5 day) midlatitude ridge over the Atlantic Ocean In early January 
1979 was energet Icullv very similar to the development of the high latitude 
Rex blocks In the same region that preceded and followed it. The similarity 
appeared in both the spectral and limited area energetics calculations. 

In all three developments, marked baroclinic energy conversion by cyclone- 
scale waves over eastern North America snd the western Atlantic preceded 
strong nonlinear interaction between the le waves and the planetary-scale 
waves. However, In two cases Rex blocking resulted, but In the third It 
did not. The explanation of the different response of the regional flow 
field to apparently similar forcing lies In the antecedent structure of the 
regional flow (as represented by the zonal mean flow plus the planetary 
waves). In the 1978-79 examples, the Rex blocking was preceded by a weak 
large-scale ridge over the Greenland-Iceland area and the absence of a 
large-scale low over eastern Canada. The midlatitude ridge development was 
preceded by a well developed Hudson Bay low. 

Simmons and Hoskins (1978, 1980) used a system with one zonal harmonic 
wave (wavenumber 6 or 9) and a variable mean flow to show that the structure 
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of ttM atan flow can affact tha growth rate and mxlnuM aaplltude of baro- 
clinic wavaa through barotroplc of facta. The Initial growth of the wavea 
In their node! la barocllnlc, but aa tha wavaa occlude tholr barotroplc 
dacay rate la profoundly affected by tha natura of tha naan flow In which 
they are enbedded and with which they Interact. 

In the real atnoaphere, the growth and decay of barocllnlc wavea are 
probably alao profoundly affected by barotroplc Interactlona with the mean 
flow within a given aector. The atructure of the flow field given by the 
axially aymmetrlc flow and existing planetary wavea detemlnea the growth 
rate, amplitude, and decay rate of barocllnlcally unatable wavea through 
modulating barotroplc interactlona. These Interactions can uc with the 
zonal mean flow or with the planetary waves. The spectral results of Kao 
and Chi (1978) and Tsay and Kao (1978) show that nonlinear barotropic 
processes play an Important role In determining the growth and decay 
rates of both large-scale and cyclone-scale waves. 

Our study shows that strongly amplified cyclone waves can have a pro- 
found effect on the development of ultralong waves through nonlinear 
interaction. The nature of the regional flow field in which finite ampli- 
tude cyclone waves occur not only has an effect on how the cyclones 
develop, but also on how they ultimately alter the regional flow field 
through interactions with it. This interaction Involves both barotropic 
energy transfer and significant sensible heat fluxes. The result can be 
the development of a short duration mid-ocean ridge with accelerated 
westerlies across the North Atlantic, or development of a blocking high 
over Greenland and Iceland that persists for nearly two weeks. 
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The growth of the blocking high over Alaska was accompanied by 
barocllnlc development of an ultralong wave trough ever East Asia and the 
western Pacific. The western Pacific large-scale trough develop- 
ment occurred over a region of large, mean SST gradients although the 
magnitude of the SST gradient In the North Pacific Is typically about one 
half of that In the western North Atlantic (Sanders and Gyakum, 1980). 
Siberian snow cover can further enhance the atmospheric barocllnlclty In 
this region. This would lead to more vigorous cyclogenesis In this region, 
Just as It does over the Atlantic (Sanders and Gyakum, 1980). Although It 
should be checked more carefully. It appears from our analysis that the 
zonal scale of the major barocllnlc disturbances In the middle and upper 
troposphere Is larger over the Pacific than over the Atlantic. Thus, the 
major cyclonic developments over the Pacific are described primarily by 
the ultralong waves, whereas Atlantic cyclogenesis Is described by 
intermediate-scale waves. Enhanced barocllnlclty In either region due to 
SST anomalies and/or high surface albedo due to unusually heavy snow cover 
Increases the activity of the dominant energetic scale In that region. Thus 
Pacific blocking can result from the barocllnlc amplification of ultra- 
long waves that Is stimulated by an enhanced barocllnlc zone in the western 
Pacific. 

Of course, large-scale barocllnlc development over eastern Siberia 
and the western Pacific does not necessarily lead to the development of 
Alaskan blocking. In fact. It can accompany the demise of Alaskan blocking 
as In the January 1979 case. Some additional factor Is needed to allow 
the large-scale barocllnlc energy conversion to create a blocking ridge. 
This may be related to the. interaction of the finite amplitude barocllnlc 
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wave idth the wean flow in wiilch It la aWbedded juet aa in the Atlantic 
caee. The 1978-79 Pacific blocking waa apparently triggered by retrogrea- 
alon of the blocking high acroaa northern Canada. However t it ia not 
clear that thla particular behavior ia a neceaaary condition for Pacific 
blocking. Other factora, auch aa the role of the atationary forcing over 
the Pacific, ahould alao be considered. 

We would expect that large-acale atationary forcing ia a more dominant 
factor over the Pacific Basin than other the Atlantic. This is because 
the Pacific Ocean has a much larger zonal extent and a much larger surface 
area than the Atlantic. In addition, the Himalayan Massif provides a 
larger upstream orographic obstruction to the mean flow than does the 
North American Cordillera. As a result, blocking theories that Invoke the 
stationary forcing of orography and land-sea differential heating as the 
mechanism for block development should be more applicable to the Pacific 
region. 

Tung and Llndzen (1979) suggested that blocking could result from the 
resonant amplification of a free Rossby wave that becomes stationary at 
a position where it is in phase with the stationary forcing of orography 
and differential heating. This proposed behavior Is at least superficially 
similar to the 1978-79 Pacific case, where t »e blocking ridge retrograded 
until it was over Alaska and in phase with the eastern Pacific standing 
wave ridge. It then amplified rapidly due to barocllnlc processes In con- 
junction with the deepening western Pacific trough. However, the growth 
rate of this Alaskan ridge may have been too rapid to be explained by 
linear resonance. Many more case studies would have to be examined to 
determine if the 1978-79 behavior Is typical of Pacific blocking. Tung 
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and Llndzen's theory may be more favorably applied to longer period circu- 
lation featurea auch aa the very peraiatent eaatern Pacific ridge in the 
winter of 1976-77 (Chen and Shukla, 1981). 

Interannual varlationa in the stationary forcing caused by alterations 
in sea surface temperatures can affect the intensity of baroclinic cyclo- 
genesis and thus the strength of the mechanisms that force block develop- 
ment (by changing the SST gradients), as well as the structure of the 
standing waves. Mean conditions promoting explor/lve cyclogenesis should 
also favor more frequent and stronger blocking. The observed preponderance 
of blocked weather patterns during certain winter seasons as opposed to 
others may be explained as a result of these variations. 

7. Conclusion 

From our spectral energetics analysis of winter season blocking case 
studies, we have found two different mechanisms through which baroclinic 
processes can lead to the development of blocking: 1) the nonlinear forc- 

ing of ultralong waves by intense, baroclinic cyclone-scale waves, and 
2) baroclinic amplification of planetary-scale waves. The second mecha- 
nism can be augmented by barotroplc transfer of mean flow kinetic energy 
to the planetary-scale eddies. The energy source for the blocking comes 
from the local region or just upstream of the ridge. The intermediate- 
scale cyclogenesis that accompanied the development of the Greenland block 
in December 1978 may have been fostered by abnormally strong SST gradients 
near Newfoundland. An enhanced baroclinic zone over the western Pacific 
may also have stimulated the large-scale cyclogenesis that accompanied the 
Alaskan block’s development, but this is less certain. 


3> 


In addition. It appears that at least two necessary conditions nust 
be Mt to allow a block to forw. One of the above forcing wechnnlsas with 
a large source of potential energy (Ajj) Must be In operation, and the 
orientation of existing planetary-scale waves Must be such that the forcing 
mschanlsw can amplify them into a blocking pattern. 

Finally, we see two major Implications of our results to numerical 
weather prediction and extended range forecasting. We have found that Che 
forcing of the ultralong waves by cyclone-scale waves can be very strong. 
Furthermore, the nonlinear interaction between Intermediate and large-scale 
waves in the atmosphere is not continuous. It occurs in discrete episodes 
when the rate of energy transfer Is very large. This Is in contrast to the 
baroclinic energy source for the ultralong waves which maintains a steadier 
rate. Thus, errors in the prediction of cyclone-scale waves will lead to 
errors in the forecasting of ultralong wave development (Gall et al., 1979). 
Models must be able to predict explosive cyclogenesis If they are to be 
able to properly forecast the development of the ultralong waves and block- 
ing. In addition, the models must accurately handle the position and phase 
of the antecedent planetary-scale waves so that they can correctly forecast 
the effects of Intense nonlinear or baroclinic forcing on the evolution of 
the existing long wave pattern. These are areas where existing forecast 
models Itave exhibited a marked deficiency. 
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APPENDIX 


The complete formulation of the apectral energetics equations used In this 
study are given by Hansen (1981). The only difference between our formulation 
and that given by Saltsman (1970) is In the nonlinear terms. Saltsman's 
flux term should be multiplied by 1/2* and the remaining half of the original 
term should be added to the nonlinear interaction term (Kanamitau et al., 1972) 
The resulting formulation of the nonlinear interaction terms are given below. 


In the spectral kinetic energy equation 
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rtprcecnta fcha rata of changa of imvanuabar n klnatlc anargy diia to nonllnaar 
tranafar of klnatlc anargy to wavamaRbar n from all other navaa. The aonal 
hamonlc axpanalon van truncated at wavanuidiar N. U» V. and n refer to the 
Fourier coafficlanta of u» v. and w, d ia the lat'*'.ude, p la the praaaure, 
and a ia the radiua of the earth. Whan aunmad over all wavenuiabarB, 1 to N» 
thla fociiulation of C|^(n|ai» n4«) yialda no net change in tlw total eddy 
kinetic energy. 
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ia the nonlinear transfer of APE from all other waves to wavenumber n. o is 
the static stability* R is the gas constant* is the specific heat of air 
at constant pressure* and T is the Fourier coefficient of the temperature. 
Note that the contribution to C^(n|m* n^) given by the fourth term 

T fi jj + T n 1 

C P ml -n n-m n -n-mj 
- P 

does not sum to zero when summed over all wavenumbers whereas contributions 
from all the other terms do sum to zero. 
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Figure Captions 

500 Kb heighta and teKperaturea at 0000 GMT 21 December 1978. 

500 mb heightc and temperatures at 1200 GMT 25 December 1978. 

500 mb heights and temperatures at 1200 GMT 29 December 1978. 

Spectrally filtered Hovmoller diagram of the planetary-scale height 
field et 70*N for December 1978 (contour Interval Is 100 m, 

positive ^values are contoured with solid lines, negative values 
are contoured with dashed lines). 


Temporal evolution of the vertically Integrated (1000-100 mb) 
energy cycle for 30*N to 80*N for 14 to 31 December 1978. 

Vertically integrated energy conversions in the same domain dur- 
ing late December 1978, 

Planetary-scale (wavenumbers 1-4) energetics for 14-31 December 
1978. The upper part of the figure shows the evolution of K, . 
(solid line) and A._. (dashed line) and the lower part shows^' 
C(A„, A._.) (heavy^solld line); C(A,_., Kj_.) (dashed line); 

C(K 2 t (dot-dashed line) and tne^contributlon to K. . due 

to nonlinear interactions of waves 1-4 with all other waves (thin 
solid line). 


Intermediate-scale (wavenumbers 5-10) energetics for 14-31 
December 1978. The upper part of tlie figure shows 
line) and Ac_|n (dashed line), and the lower part shows C(A^, 
A-^ij,) (thicK solid line); C(A^ ,^) (dashed line); 


C(Kj 5 » (dot-dashed line);^and the rate of change of K 
due 'to nonlinear interactions with all other waves (thin so 
line) . 


!iP 


The intermediate-scale height field (Zc.ia) at 500 mb at 1200 
GMT 18 December 1978 (contour interval^ls 100 m, solid contours 
indicate positive values and dashed contours indicate negative 
values in this and all other polar sterographic projections). 



at 500 mb on 1200 
5 X 10“3 m^/s^. 


GMT 18 December 


1978 


(contour 


The intermediate-scale heat transport (vT, ._) at 850 mb on 
1200 GMT 10 December 1978 (contour interval'*^ls 50*C m/s). 

The intermediate-scale momentum transport (nvc_,Q) at 300 mb on 
1200 GltT 19 December 1978 (contour interval la 100 m^.'s^) . 

at 500 mb on 1200 GMT 19 December 1978. 
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Figure 10a. Kc.ia at 500 mb on 1200 GMT 19 December 1978 (contour Interval 
la l60 m2/a2). 

Figure 10b. K, . at 500 mb on 1200 GMT 23 December 1978 (contour interval 
ai in a) . 


Figure 11. 
Figure 12a. 
Figure 12b. 
Figure 13. 

Figure 14. 


C(A^_^, K^_^) at 500 mb on 1200 GMT 19 December 1978. 

*1-4 ^200 GMT 23 December 1978 (contoura every 100 m). 

at 500 mb on 1200 GMT 23 December 1978 (contoura every 5*C). 

vTi «4 ^200 GMT 21 December 1978 (contours every 50*C 

m/s) . 

Limited area energetics for the region bounded by 85*W, 20*W« 

30*N and 75*N. 


Figure 15. 
Figure 16. 
Figure 17. 
F( ' ire 18. 

Figure 19a. 

Figure 19b. 
Figure 20a. 

Figure 20b. 

Figure 21a. 
Figure 21b. 
Figure 21c. 
Figure 21d. 
Figure 22. 
Figure 23. 
Figure 24. 
Figure 25a. 
Figure 25b. 


500 mb heights and temperatures at 1200 GMT 30 December 1978. 

500 mb heights and temperatures at 0000 GMT 4 January 1979. 

500 mb heights and temperatures at 0000 GMT 12 January 1979. 

Spectrally filtered Hovmoller diagram of . at 60*N for late 
December 1978 to mid January 1979. 

Energy cycle for late December 1978 to mid January 1979 (other- 
wise same as Figure 5a). 

C(A^t Ag), C(Ag, Kg) and C(K2, Kg) for the same period. 

Planetary-scale energetics for 24 December 1978 to 12 January 
1979 (otherwise same as Figure 6a). 

Intermediate-scale energetics for the same period (same as 
Figure 6b otherwise). 

*1_4 at 500 mb on 1200 GMT 30 December 1978. 
vT^_^ at 500 rob on 1200 GWT 30 December 1978. 

T^_^ at 500 mb on 1200 GMT 30 December 1978. 

G(A^_^, Kj^_^) at 500 mt on 1200 GMT 30 December 1978. 

“'^X_4 ^200 GMT 30 December 1978. 

Kj_^ at 500 mb on 1200 GMT 31 December 1978. 
vTj_^ at 500 mb on 0000 GMT 4 January 1979. 

Tj^ at 500 mb on 0000 GMT 4 January 1979. 

*1-4 St 500 rob on 0000 GMT 4 January 1979. 
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ABSTRACT 

The height and temperature fields at the mandatory levels during a two- 
month period (December 1977-*January 1978) are used to evaluate the conversion 
between available potential energy (A) and kinetic energy (K)« C(A,K). The 

results reveal: 

1. The energy conversion due to the mean meridional circulation, 

A 

C(A 2 « Kj,) ■ - 0.18 and due to the eddy motions, 

C(Ag, Kg) - 2.6 Wm"^. 

2. ^(Ag, Kg) is attributed to the standing (2-month mean) mode, while 
Kg) is due to the transient eddies. 

3. The contributions to C(Ag, Kg) from various wavenumber regimes, 
long waves (n*l-4), and cyclone and short waives (n*3-15), are: 

a. Transient cyclone and short waves are dominant in the midlatitudes 
and the middle layer of atmosphere. 

b. Transient long waves are significant in the middle layer of tie 
atmosphere at mid- and high latitudes. 

c. Standing long waves are not negligible in the mid latitudes and 
the middle layer of atmosphere. 
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Introduction 

The conversion bettfeen available potential energy (A) and kinetic energy 
(K), C(A, K), is regarded as the baroclinic process that maintains the atmospheric 
circulation (Ulln-Nlelsen, 1968). Tomatsu (1979) recently made an extensive 
computation of this energy conversion employing multilevel data. However, 
he has noted that there had been very few analyses of C(A, K) for the hemisphere. 
This might well be due to the difficulty of evaluating the u field in the 
atmosphere. 

Recently, we used an alternative mathematical algorithm to solve the u) 
equation (see Appendix) rather than the conventional relaxation method. Some 
interesting features of C(A, K) were found in our computations using a two- 
month (December 1977-January 1978) multilevel data set analyzed at the National 
Meteorological Center (NMC). The purpose of this note is to report these 
C(A, K) results in the wavenumber regime and the contributions to C(A, K) 
from the standing and transient modes of the atmospheric motion. 


Computation and Data 

The conversion between A and K, 


C(A, K) 


I 

8 


f _ R yxdsdp, 

Jo Js P 


can be expressed in terms of a Fourier expansion representing the contribution 
from various wave components. 
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C(A, K) 


1 

g 


fPr 


0 ■'gj 


•^0 
0 J 


C(0) 


T“ — (u)T + uT) cos^d^dp 
p nn nn rri- 

N 

I C(n) 

n"l 


( 1 ) 


C(A^. K^) 


C(Ag, Kj,) 


Till! notations uHtui are conventional and ( ) ® complex conluj’.ate ot ( ). 

The two terms on the right-hand side of (1) are the energy conversions due 
to the mean meridii tuil and eddy motions, respect Ivelv. 

The expression of c;(A, K) can also be resolved lurther to examine the 
contributions from the standing and transient modes of the atmospheric circulation, 


C(A.- K) - ‘ - cobMWp + i 

g jo P 0 0 g 


fPoi 


' p‘^0^ cosiJ>d(lidp 


C(0). 


C(0),j, 


I fPo 

g Jo i 


C(A^, 




C(A^. 


( 2 ) 


) — (u! T 4- ( 1 ) T ) cos(J>d(|)dp + 


'4> n-1 


pnn n n 


^0 


>0 
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T --(..)'T’ + ui'*T') cos^id^idp 
^.p nn nn 
<J) n*l 


I C(n), 


n~l 


N 

I C(n), 
n*l 


C(Ap^, 


C(Ag, 


where ( ) * time mean and ( )' ■ ( )-( ), is the transient mode. The subscript.*- S 
and T denote the standing and transient modes, respectively. N ■ 15 and che 
time mean is a two-month average in this study. 

The p vertical motion of the atmosphere is evaluated by using the oj equation 
of the quasi-geostrophic model. The tu equation has the form of Poisson's 
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equation and conventionally la solved by the relaxation method (e.g., 
Krishnamurti, 1968). In this study, an alternative approach employing ilie 
tech, ^que of spherical harmonics to solve the u equation is used. The. a> field 
can be obtained through a matrix equation which relates the spectral coefficients 
of u) and the forcing functions in the u equation. Our experience reveals hut 
this approach Is computationally faster than the conventional relaxation method. 
The detailed algorithm of this approach Is shown In the Appendix. The boundary 
conditions used In this study are w • 0 at 1000 and 100 mb for the simplicity 
of computations. 

The computations of this study use the twice daily height and temperature 
fields on the 2.S” x 2.3° grid of t.ie UHL analysis for December 1977 and 
January 1978. The data are at the mandatory levels; 1000, 850, 700, 500, 

400, 300, 250, 200, 150, and 100 mb. The domain of integration In (7) and 
(2) i.s from 20°N to 85°N and from 1000 mb to 100 mb. 

Re s ul t s 

An overall view of the conversion between available potential energy 

and kinetic energy for various modes of the atmospheric motion is provided in 

Table 1. Note that the eddy conversion is also divided into two groups; 

4 15 

'I C(n) represents the long-wave regime and C(n) the cyclone and short- 
n®l n*5 ^2 

wave regime. The numerical values of C(A^, Kjj) ■ - 0.18 Wm and K^) * 

2.6 Wm"^ in the present study. Tomatsu obtains a two-month (December and 

— _2 — -2 

January) mean values of C(Aj,, K^) “ - 0.6 Wm and C(Ajj, Kg) • 3.14 Wm . Our 

results are somewhat smaller than Tomatsu’ s, but not unrealistic. In fact, our 

results are also close to Oerleraans' (1980) recent estimate of C(Ag, K^,) 

“• -2 

which is 2.1 Wm in winter. The salient features shown in Table 1 are: 
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4 

(1) ) C(n) and ) C(n) are comparabje; (2) C(A^, K,^) i- uttrlhuleti to the 

n*l n*-j 

standing mode and C(A„, Kp) is mainly contributed by the transient eddies; 

^ 4 

(')) the major amount of T C(n) is due to the transient eddies, but the con- 

n»l 15 

tribution from the standing eddies Is not insignificant, and; (4) I C(n) Is 

n«5 

raulnlv provided by the transient eddies. 

The spectral d Intr Ihiit ion of C(n) is shown in Figure 1. The maximum value!! 
C(n) In this study appear at wave numbers 2 and 6. This is in agreement with 
the spectral energetics results compiled by Saltzman (1970). Tomatsu's com- 
putation shows that the maximum C(n) occurs at waves I and 3 in winter. The 
hatched area of Figure 1 represents the spectrum of C(n)^. The maximum value 
of C(n)^ occurs at wave 6. The contribution of C(n)g comes from waves 1-3. 

In order to have a better view of the climatology of the conversion 

between A and K, Figure 2 displays the latitudinal height distribution of the 

energy conversion for various modes of atmospheric motions and various wave- 
number regimes. Tht distribution of C(Aj.. K^), Figure 2a, is similar to that 
shown in Tomatsu’s study. The negative values appear between 50"N and 70“N 
which may be the result of the Ferre! circulation. The positive values occur 
between 35“N and 50*N, and in the polar are«, north of 70“N. Although the areal 

average of CfA^, K^) is smaller in this study, the numerical values in the 

C(A^, K^,) distribution in the current study are larger than Tomatsu’s. This 
indicates that the intensity of the mean meridional circulation in the winter 
analysed in this study is stronger. 

The distribution of C(Ag, K^), Figure 2b, shows three areas of maximum 
value at 500 mb: 4S“N (major one), 65“S and 75“N. Tomatsu’s result only 

shows a maximum value at 600 mb and 50"N. It is of Interest to examine the 
contribution of different wave regimes to C(Ag, K^) . The comparison between 
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4 i:> 

C(A„, K„), I C(n) (Figure 2c) and J C(n) (Figure 2d) ahowa very clearly that 
the major maximum of C(A^, Kg) at 45*N la mainly attributed to the cyclone and 
short waves. The two minor maxlM of C(Agi K^) at 65*N and 75*N are due to the 
long waves. This may Indicate that the wave activities in the high latitudes 
are dominated by the long waves, while in the mid-latitudes cyclone and short 
waves are active, in addition to the long waves. 

In the right-hand column of Figure 2, we display the latitude-height 
distributions of either the standing or transient modes of those four 
quantities discussed above depending upon which provides the dominant con- 
tribution to their totals, respectively. The C(Ag, Kg) distribution is more 

or less decided by C(Ag, ^2^8* distribution is determined 

4 

by C(Ag, Kp)'j.- The further breakdown of C(Ag, Kg).j, Into 1 C(n).j, and 
15 ' 15 15 

I C(n)-. shows that I C(n) distribution is almost identical to I C(n) . 
n*5 4 n*5 4 n*5 

However, ][ C(n) and ][ C(n)„ has some significant difference over the loca- 
n"l n*l 

Kg) occurs, l.e. 500 mb and 45*N. The 
difference between these two quantities indicates that the major contribution 
of the standing eddies appears at 500 mb and 45*N. 


The idea of computating the energy conversion between A and K is not new. 


However, very few studies have been made to evaluate this energy conversion 
using multilevel observational data. Some interesting features shown in our 
C(A, K) computations may be useful in evaluating model simulations 
of the atmospheric general circulation. 
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Table 1. Contributions to conversion between available potential 
energy (A) and kinetic energy (K) from various modes 
of atmospheric motions. Units Wm**^. 


Total (S4-T) Standing (S) Transient (T) 


K^) 

-0.18 

-0.22 

0.04 

C(Ag, Kg) 

2.60 

0.40 

2.20 

4 

1 w 

n*J 

1.15 

0.36 

0.79 

15 

I C(n) 

1.45 

0.04 

1.44 


n*5 
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Apptndix 

The u equation of the quasi-geoerrophlc nodel nay be vritten aa 


V^w 


5»^u> 

17 


f . 


where ^ “ f •*^'*'* ^ ~ J(\|.. |J-) . 


(Al) 


II ■ i)(p) and ♦ i« geopotential. The forcing of function in (Al) coniaina only th»> 
vertical differentiation of vorticity advectlon and the Laplacian of thermal advec- 
tloii. last u« expreMH ui and F In lernut of aolid spherical harmonics (I'laizman, 


& k 


in 


U) 


m 

'n 


(AJ) 


and 


F - . >; f" y"* 

n n n 


(A3) 


wliere Y™ “ p"’ e^*”'" .oui P^fu) is the Associtited I.epcndre function ot the 
n 11 n ' 

latitude, <^, at which p ■ ein^. and F™ are the spectral coefficients 
of ( 1 ) and the forcing function, respectively. Substitution (A2) and (A3) 
into (A4), we can obtain the spectral form of (A2), 


n(n-4j) 
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whero a ia the earth's radius. The second-order derivative of u" in (AA) 

n 

is evaluated by the finite difference scheme illustrated in Figure Al. 

The second-order derivative of w® at the level K is evaluated by dividing 

the difference between the first-order derivative of w"* in the middle of 

n 

levels (K + 1 , K) and of levels (K, K-1) by ^ Therefore, 

(AA) can be written as, 


Vl % + X) - a,™ (K) + (K-1) - 


f” (K) 


(A5) 


where 


2f' 


'K+1 


.2 


n(n+l) 


2f‘ 


0 Apj^ 


2f 


'K-1 0 APj^ ^^Pr+1 *Pr^ * 


Applying (A5) to every level, we can obtain a linear system which can be 
expressed in terms of a matrix equation. 


where 
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M u * F 
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• • 

• • 
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N l8 tha flMxlMUit lUMbar of vortical lovela, R la a trladlagonal aquare matrix 
and to)" at the varloua levela can be obtained by aolvlng 

The procedurea are aa followa: (a) calculate the forcing function of 

the to) equation on the 2.S* x 2.5* grid; (b) uae Ellaaeaaer'a (1966) acheme 
to create the spectral coefflclenta of the forcing functions; (c) the spectral 
coefficients of the vertical motion are solved by equation (A6); and (d) the 
to) field is constructed by equation (A2). 

The spectral truncation used In this study Is triangular with 0 |m| < 25 

and 0 < n < 25. Since only hemispherical data Is available, we roust make 
some assumption concerning the parity at the physical variables. We assume 
that to) and P are even. 


■iirrs— 
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Figure 1. Spectral dtetrlbutiona of C(n) (solid line), €(n)^ (solid line 

enclosing the hatched area) and C(n)^ (dashed line). I'nlt: 

-1 >2 
10 * Wn 

Figure 2. Latitudinal-height distribution of (a) C(k^, K 2 ). (b) C(A^, K^), 

(c) I C(n), (d) I C(n), (e) C(A-, ILL. (f) C(A-, 0„. 
n-l n»5 

4 15 

(8) X C(n)^ and (h) X C(n)^. Unit? lO"^ Wm ^ mb’*. 
n"l n«»5 

Figure A-1. Schematic diagram for the computation of the vertical derivative 

for u* . 
n 
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Abatract 

In January 1977* a waning of tha polar tropoaphera occurrad following 
a major stratoapharlc wanning. Tha tomparatura incraaaa in tha polar tropo- 
sphere ia found to be larger below 400 mb than in the upper troposphere. 
Analysis of the heat h’idget of the polar troposphere shows that the eddy 
heat flux convergence ia the major factor responsible for tha tropospheric 
warming. 

Introduction 

During January 1977, a major sudden warming in the stratosphere was 
followed by a striking temperature Increase in the polar troposphere (Quiroz, 
1977). This sequence of events raises the question of causal connection 
between these two occurrances. It would be interesting to investigate whether 
there is a direct connection between the stratospheric and tropospheric warm- 
ings. Recently, Koermer and Kao (1980) studied the evolution of the kinetic 
and thermal energy as well as the heat and momentum flux convergence in the 
stratosphere and troposphere during this warming compared to the minor warming 
of 1975-76. O’delll and Taylor (1979) investigated the 20 mb and 300 mb heat 
and momentum budgets and found chat warming in the polar upper troposphere in 
early January was due to increased subsidence that was forced by southward eddy 


momentum transport. 
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We have computed the heat budget in the tropoaphere %rhere we found that the 
largeat temperature increaee over the polar cap in mid-January 1977 occurred below 
400 mb. Taylor and Perry (1977) have pointed out that an unuaually 
large number of 1000-500 mb thickneaa anomoliea centered around 80*N latitude 
occurred during mid-January of thia year. The purpoae of thia note ie to show 
that the lower tropoapheric warming was primarily the result of strong eddy heat 
flux convergence in the lower troposphere. Although a complete explanation of 
thu cause of the enhanced oddy transport of sensible heat into the polar 
tropoaphere is beyond the scope of the present report* a possible mechanism 
Is proposed. 


Analysis Method and Data 

The heat budget of the polar trophosphere is investigated by using the 
zonal ly averaged form of the thermodynamic equation: 



1 X 

a cos 3(^ 






’!!£ . -R. T 






9p c p 

w - 

z 

- 

3P ~ Cpp 



The z subscript denotes variables that wero averaged around a latitude circle 
and the E subscript denotes eddy quantities that were departures from the 
zonal average. The term on the left-hand side is the local rate of change of 
the zonally averaged temperature. The terms on the right-hand side are the 
temperature change due to 1) eddy heat flux convergence: 2) convergence of the 
zonal mean heat flux: 3) adiabatic warming (or cooling) due to zonal mean 
vertical motion; 4) adiabatic warming due to eddy vertical motion; and 
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S) dlabatic heating. 

Data friM the ubjective analyala achana of the National Meteorological 
Center (NMC) for ovary 2.5* of longitude at the 10 Mandatory lavala (1000, SSO, 

700, 500, 400, 300, 250, 200. 150, and 100 nb) in the tropoaphera are uaed to 
evaluate the varioua tame in tnia aquation at 70*N and 80*N. NW: enployed the 
Flattery analyaia achena (Flattery, 1970) in 1976^77 which renovea the divergent 
component of the wind from the horiaontal notion field. ■ 0 in the Flattery 
analyaia ao the raaidual in the heat budget repreaenta not only the effecta of 
diabatic heating but alao changea in T 2 due to the actual convergence of the zonal 
mean heat flux ax well aa any effecta due to errora in the data. The exact quality 
of the d^^a at polar latitudea ia difficult to detemine, but for the aake of 
dlacuaaion, chatigaa in or contrlbutlona to changea In with nagnitudea of l*K/day 
or leaa nay be conaldered to be noiae. The vertical notion ia conputed from the 
quaai-geoatrophic form of the u-equation aaaunlTig 0 ■ 0 at 1(X)0 nb and 100 mb 
(Chen et al., 1981). The unequatlon includea the effecta of both dynamic and 
thermodynamic forcing on the vertical notion field (Wiin-Nlelaen, 1964). Theae 
include the vertical differentiation of the abaolute vorticity advection and the 
Laplaclan of the themal advection. The oi valuea obtained from thia method may 
be undereat Imated to aome extent due to exceaaive amoothing of the height field 
by the objective analyaia acheme in data aparae regiona (Wiin-Nielaen, 1968). 

Dlacuaaion and Speculation 

Figure I ahowa the average croaa-aection of the zonal mean temperature for 
10>18 January 1977. The following dlacuaaion will concentrate on the reaulta 
at 80*N becauae the largeat tropoapherlc temperature riae occurred at or north 
of 80*N (Figure 1). 
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The xonal meim tenperacure ac 80“N aa a function of haighc and time la 
ahown in Figure 2. The largest; tenperaCure rise occurs In the lower tropo- 
aphev'e, below 700 mb, from 9 to 12 January with a slightly smaller tempera- 
ture Increase in the mid-troposphere. Very little temperature rise la 
i>vident In the upper troposphere (300-200 mb). The largest rate of temperature 
increase occurs at low and mid levels (Figure 3a), but warming extended 
throughout the troposphere below 300 mb from 1200 GMT 9 January through 
0000 GMT 12 January. Significant heating due to eddy heat flux convergence 
occurs In the layer below 400 mb from 1200 OfT 9 January through 0000 GMT 
12 January (Figure 3b) and colncldea almost exactly with the observed tem- 
perature Increase. The maximum warming due to eddy sensible heat transport Is 
4*K/day between 700 and 1000 mb. Comparison of Figures 3a and b shows 
that the eddy heat flux accounts for almost the entire observed temperature 
rise in mid-January. A similar relationship between the rate of change of 
the zonal mean temperature and positive eddy heat flux convergence In the 
lower troposphere is also apparent In the heat budget results at 70*N 
(not shown). 

Warming due to zonal mean adiabatic descent at 80*N Is present through- 
out a deep layer on 8-9 January with a maximum value of 3*K/day between 200 
mb and 300 mb (Figure 3c). However this effect Is not reflected In a rise 
in T^. The adiabatic warming in our calculation la compensated by the 
residual (Figure 3d). The effects of vertical motion are unimportant at 70“N 
(not shown). The large eddy heat flux convergence In the lower troposphere from 
9 to 12 January would force upward vertical motion and diabetic cooling. How- 
ever, subsidence Is indicated In Figure 3c on 9 and 10 January. Apparently, 
the strong southward eddy momenttun flux occurring at this time at polar latitudes 
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(O'Neill and Taylor« 1979) forced deaeending notion that dominated the effect of 
the northward eddy heat flux on the mean vertical velocity. The ineignif icance 
of adiabatic warming together with the absence of a large temperature rise In 
the upper troposphere (Figure 2) indicate, as O'Neill and Taylor (1979) have 
already shown, that warming due to adiabatic descent did not propagate downward 
from the stratosphere into the troposphere. Eddy vertical motion did not 
significantly effect the polar heat budget. 

Koemer and Kao (1980) note that once the 1976-77 major warming was in 
progress, there was a decoupling of the stratospheric and tropospheric circula- 
tion and blocking in the troposphere. Under normal conditions, the band of 
strongest zonal mean westerly winds arts ns a waveguide for vertically propagating 
planetary-scale waves (Dickinson, 1968). During a major stratospheric warming 
the westerly jet in the stratosphere Is replaced by mean easterlies preventing 
the normal upward propagation of planetary waves (Charney and Drazin, 1961). 
Although easterly winds appeared almost simultaneously at polar latitudes in 
both the lower stratosphere and troposphere shortly before 10 January 1977 
(Quiroz, 1977; O'Neill and Taylor, l'*79), substantial warming in the middle 
and upper stratosphere occurred from 28 December 1976 to 3 January 1977 . 

(Quiroz, 1977). Reversal of the zonal wind and temperature gradient occurred 
in the upper stratosphere before extending downward (Koermer and Kao, 1980). 

This warming pulse in late December and early January may have inif.i <ted tho 
disruption of the westerly wave guide in the stratosphere. As a r. lult, ampli- 
fying tropospheric planetary waves associated with Atlantic blocking in the 
first 10 days of January may have been unable to propagate vertically. The 
tropospheric planetary waves reached their maximum amplitude at about the same 
time or sho«*tly after easterly winds appeared in the lower stratosphere (Table 1) . 
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The easterly mean flov In the lower stratosphere would completely disrupt 
the westerly wave-guide. The coincidence of tropospheric wave amplifica- 
tion with the appearance of stratospheric easterlies may have lead to the 
observed eddy heat flux convergence at polar latitudes. With the waves 
unable to propagate upward away from their source region, the northward heat 
transport normally associated with vertically propagating waves may have 
manifested Itself in th# large observed sensible heat transport Into the polar 
troposphere. The heat transport accompanying the amplifying waves lead to 
warming near the pole Instead of upward wave propagation. 

McGulrk (1977) noted that major perturbations of tha zonal mean tropo- 
spheric temperatures with high latitude warming and mid-latitude cooling 
occur after the commencement of stratospheric warmings. Labltzke (1965) 
has shown that tropospheric blocking patterns occur roughly 10 days after 
the onset of stratospheric warmings. I blocking, a pronounced warming 
occurs In the northern part of the blocked zone with cooling in the southern 
part (Rex, 1950). Consequently, a mechanism similar to one presently pro- 
posed may have been responsible for the growth of the blocks studied by 
1 /ibltzke and for the pattern of T_ changes noted by McGulrk. 

mt 

In conclusion, the tropospheric polar warming In January 1977 was 
principally the result of eddy heat flux convergence In the middle and 
lower troposphere. It appears that the tropospheric warming was not due to 
adiabatic descent within the troposphere. This warming may have been a 
manifestation of the crapping of vertically propagating planetary waves 
beneath stratospheric easterlies. Thus, the stratospheric warming may have 
had an indirect effect on the tropospheric circulation. 


Acknowltdgg— nt ■ 


This study vss supportsd in psrt by tha Nationsl Asronautics and Space 
Administration under Grant NSG-S339, and the National Science Foundation 
under Grant ATM-7924568 • 


T aaMwausiSBBg S S . 




Table 1. Dates of the maximum 500 a amplitude of the sum of 
zonal harmonic height waves 1 through 4 over the 
Atlantic In early January 1977 


Latitude 

Date 

Maximum 
Wave Amplitude 

Longitude 

50*N 

8 Jan 

475 m 

25 W 

60*N 

iO Jan 

419 m 

30 W 

70*N 

11 Jan 

378 m 

40 W 


FIGURE CAPTIONS 


Figure 1. The average meridional croaa eectlon of T^(*K) for 0000 GMT 10 
.T luary 1977 through 1200 GMT 18 January 1977. 

Figure 2. Ihe ..onally averaged temperature, T 2 (*K). at 80*N aa a function 
of height and time. 

Figure 3. Height-time aectlona of the terma in the heat budget at 80*N: 

a) the time rate of change of the sonally averaged teimierature, 

1 3 

3Tjj/3ts b) eddy heat flux convergence, - ^ 

c) adiabatic varmlng (cooling) due to xonal mean vertical motion, 
and; d) the reaidual. The units are *K/day with positive values 



shaded . 


REFERENCES 


Chon, T.C., A.R. Hansen, and J.J. Trlb*^ia, A note on the release of available 
potential energy. Submitted to J. Meteorol. Soc. of Japan, i981. 

Dickinson, R.E., Planetary Roaeby waves propagating thtuugh ;^4ak westerly wind 
wave guides. J. Atmos. Scl., 25, 269-279, 1968. 

Flattery, T., Spectral models for global analysis and forecasting. Proceedings 
of the Sixth AWS Technical Exchange Conference, U.S. Naval Academy, Air 
Weather Service Technical Report 242, pp. 42-54, 1970. 

Koermer, J.P., and S.K. Kao, Major and minor stratospheric warmings and their 
Interactions on the tropospb'^re . Pageoph, 118, 428-451, 1980. 

Labltske, K., On the mutual relation between stratosphere and troposphere 
during periods of stratospheric warmings in winter. J. Appl. Meteorol. 

4, 91-99, 1965. 

McOuirk, J.P., Planetary-scale forcing of the Janaury 1977 weather. Science 
199. 293-295. 1978. 

O'Neill, A.. f»nd B.F. Taylor. A study of the major stratospheric warming of 
1976/77. Quart. J. Roy. Meteorol. Soc., 105, 71-92, 1979. 

Quiroz, R.S., The tropospheric-stratospheric polar vortex bveakdown of January 
1977. Geophys. Res. Lett., 4, 151-154, 1977. 

Rex, D.F. , Blocking action In the middle troposphere and its effect on regional 
cHmate. I. An serological study of blocking action. Tellus, 2, 196-211, 
19S0. 

Taylor, B.F. , and J.D. Perry, The major stratospheric warming of 1976-1977. 
Nature, 267, 417-418, 1977. 

Wlln-Nlelsen, A. , On the Intensity of the general circulation of the atmo- 
sphere. Rev. Geophys., 6, 559-379, 1968. 

Wlln-Nlelsen, A., On energy conversion calculations. Mon. Wea. Rev.. 92, 161-167. 
1964. 



Latitude 





















